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SELECTION OF OPTIMAL SITE FOR NEW DEPOT OF 
SPECIFIED CAPACITY WITH TWO OBJECTIVES 


SATYA PRAKASH 


Department of Mathematics, Birla Institute of Technology & Science, Pilani 
(Rajasthan) 


AND 


ViveEK SAINI 


Software Engineer, Tractor Division, Escorts, Faridabad (Haryana) 


(Received 15 October 1987) 


The problem of selecting an optimal site for a new depot of specified 
capacity from several potential sites along with determining an optimal 
schedule to take buses from depots to the starting points of their routes and 
also determining the spare capacity available at each of the depots after the 
construction of anew depot is considered. Capacities of the respective 
depots and the number of buses required at the starting points of the respec- 
tive routes are specified. The problem has two objectives—one primary and 
another secondary. The primary objective is to minimize the capital 
expenditure to be incurred in constructing a new depot plus the present value 
of the expenditure to be incurred in total dead mileage over a planning hori- 
zon. The secondary objective is to minimize the maximum of the dead 
mileage of individual buses. An algorithm is developed to obtain the solution 
of this two-objective location of problem. 


1. INTRODUCTION 


The problem of selection of an optimal site for a new depot from several potential 
sites is of current interest. For depots built in the past cannot provide parking facili- 
ties for all buses because of increase in their number with the passage of time and 
depots have reached upper limit of augmentation of their capacities. In the past, this 
problem has been tackled by applying certain rudimentary techniques based on com- 
monsense and experience, But this has resulted at times into expenses which are 
avoidable. Owing to this and awareness that even a small percentage of savings in 
urban transportation schemes can result into substantial savings because of huge in- 
vestments involved in them, there is need to minimize unnecessary expenditure by all 
means. One of the ways in which unnecessary expenditure can be minimized is to em- 
ploy analytical tools to deal with the problems of urban bus transportation. Recently 
Sharma and Prakash? have applied analytical tools for optimizing dead mileage in 
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urban bus routes. The present paper deals with the problem of selecting an optimal 
site for a new depot of specified capacity from several potential sites along with deter- 
mining an optimal schedule to take buses from depots to the starting points of their 
routes and also determining the spare capacity available at each of the depots after the 
construction of a new depot. Capacities of the respective depots and the number of 
buses required at the starting points of the respective routes are specified. The problem 
has two objectives—one primary and another secondary. The primary objective is to 
minimize the capital expenditure to be incurred in constructing a new depot plus the 
present value of the expenditure to be incurred in total dead mileage over a planning 
horizon. The secondary objective is to minimize the maximum of the dead mileage 
of individual buses. The term dead mileage rcfers to the distance traversed by a vehicle 
when no service is provided. For instance, the distance traversed by a bus from a 
depot where it is parked over-night to the starting point of its route in the morning is 
dead mileage because this distance is traversed without severing any useful purpose. 
An algorithm is developed to obtain the solution of the two-objective location problem 
and is illustrated through a numerical example. 


2. FORMULATION OF THE PROBLEM 


Suppose that there are m existing depots, s potential sites for a new depot and 
n starting points of routes. The total number of buses required at the starting points 
of the various routes is greater than the number of buses to be parked overnight at the 
depots necessitating the construction of a new depot. Capacity of the new depot is so 
specified that there is space at depots for parking a certain number of additional buses 
besides providing parking facilities for all buses operating on routes after its construc- 
tion. For the convenience of notation, the proposed new depot at potential site / 
(1 = 1, ...,.5) will be designated as depot m+/. Let ai (i = 1, ...,m+-s) be the 
maximum number of buses that can be parked overnight at depot i, b; (j= ly oom) 
be the number of buses required at the starting point of route j, bn, be the number 
of additional buses apart from buses Operating on routes which can be parked over- 
night at depots after the construction of a new depot, cm,z (J = /, ..., s) units be the 
cost of creating parking facility for one bus at site l, day Gra 15, rn Chey SL) ke 
.e, 1) units be the distance from depot i to the starting point of route J, kms (1 = 1, 
--+, 5) units be the initial setup cost of constructing depot at site /, N be the number of 
years over which the planning horizon is spread, p units be the expenditure incurred 
in traversing one unit of dead mileage by a bus, r be the rate of interest per annum, 
m1 (1 = 1, ..., 5) be an integer assuming value 0 or | according as site / is not sele- 
cted or selected for new depot, xij (i = emg b= 5s; J = 1, ..., n) be the number of 
buses to be sent from depot i to the starting point of route J, and xi(n.1) (i = |] 
m +s) be the number of buses apart from buses operating on routes which ar be 
parked overnight at depot i after the construction of a new depot. The capacity of 


each proposed new depot at potential site is identical and is equal to "S bj — z at. It 
j>1 fa 4" 
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is required to select an optimal site for a new depot from the potential sites along with 
to determine an optimal schedule to take buses from depots to the starting points of 
their routes and also to determine the spare capacity available at each of the depots 
after the construction of a new depot subject to the constraints of the problem. The 
primary objective is to minimize the capital expenditure to be incurred in constructing 
a new depot plus the present value of the expenditure to be incurred in total dead 
mileage over the entire planning horizon. The secondary objective is to minimize the 
maximum distance among the distances traversed by individual buses from depots to 
the starting points of their respective routes. The mathematical formulation of this 
two-objective location problem is as follows: Find Amy: = Oor 1 (/ = 1,...,5) and 
xj > OCG =1, ...,.m+ s;7 = 1, ...,n + 1) which minimize 


if $ Nel ™m n } 
1S Amgz [kmgr + (DS by — EF at) emir + LF cm,yj x(m,di] 
| i=1 j=1 i=1 j=1 
C= 4 \ 
| m n | 
| +- Zz Cty X8j 
{ i=1 j=l J 
(1) 
where 
v . 
cy = 2 (365) p dig SS (1 + 7/100)-' G§ = 1, ...,m + 5;7 = 1, ..., n) 
t=1 
éee(2) 
and 
D = max {dij : xi > O(i = 1, ... m + 5;j7 = 1, ..., m)} ets) 
according to priorties in the order of their occurrence subject to the constraints 
x= Ami = | -..(4) 
l=1 
n+l 
= xij =a (i = 1,..., m) ...(5) 
j=1 
8 1 Nel m 
SE Amix(m.pi = % b- E a ...(6) 
l=1 j= 1 j=1 i=1 
m &§ : 
DS xp t SS Amyt x(myns = bj (fj = 1, ...,2 + 1). ae eh 
i=1 I=1 


3. SOLUTION PROCEDURE 


The two-objective location problem formulated above is a mixed integer non- 
linear problem. A procedure is outlined to obtain the solution of this problem. As 
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Am.,1’s are integers assuming values 0 or ! and are subjected to the constraint (4), it 
follows that all elements except one in each distinct set of values assumed by Amst, 
...» Am,s are zero and that the nonzero element is 1; so the total number of distinct 
sets of values assumed by Am4i, .. , Amys is only s. Setting Am,1 = 1 and Am,4i = 0 
(I = 2, ..., 5), the two-objective location problem given by equations (1) through (7) 
yields the following problem. Find xij > 0 (7 = lu. mm ye hy 1 ya 1) 
which minimize 


i s 0 5 Cig Xtj (8) 
Ci = kml + ( oh bj = re ai )/ Cm41 rm be 
where 
N . 
ctj = 2 (365) pay 2 (1 + 7/100) * (i = 1, ...,.m,m+1;j = 1, <eig FO 
t=1 
. (9) 
and 
Dy = max {dyy: x15 > O (i = 1, ...,.m, m+ 1; j = 1, ..., a} .. (10) 


according to priorities in the order of their occurence subject to the constraints (5), 


N+1 n+l m 

x= xmns = SF by — SF as eetiae 
j-1 jal i=1 
m,m41 

= xe = b3(j = l,....n + J). atlas 
int 


The problem described by equations (8), (9), (10), (5), (11), (12) is designated 
as the Ist two-objective problem of the given two-objective location problem. The 
problem obtained from the given two-objective location problem by setting Amy = 1 
and Am, = 0 (/ = 1, 3, .... s) is designated as the two-objective problem of the given 
two-objective location problem, and so on. Finally, the problem obtained from the 
given two-objective location problem by setting Amys = 1 and Am,1 = O (J=1,... ,s—1) 
is designated as the sth two-objective problem of the given two-objective location pro- 
blem. These two-objective problems of the given two-objective location problem 
are of similar type and are reduced to equivalent single objective transportation-type 
problems following the procedure developed by Prakash2. The single-objective tran- 
sportation-type problem equivalent to the /th two-objective problem of the given 
two-objective location problem is designated as the /-th single objective transportation- 
type problem of the given two-objective location problem. Among all the single- 
objective transportation-type problems of the given two-objective location problem, 
the one whose objective function has minimum value provides complete solution of the 


given two-objective location problem and the optimal site for a new depot is the site 
corresponding to it, 
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The procedure to reduce the Ist two-objective problem of the given two-objective 
location problem to the Ist single-objective transportation-type problem of the given 
two-objective location problem is briefly described below. First, the set {dig : i= 1, 
..m,m + 1;j = 1, .... n} is partitioned into subsets Liz (k = 1,...,q) in the following 
way. Each of the subsets Lx consists of the dij’s having the same numerical value, Li 
consists of the dij’s having the largest numerical value, Z2 consists of the dij’s having 
the next largest numerical value. and soon Finally, La consists of the dij’s having 
the smallest numerical vaue. After this, priority factors. Mo, M1, ..., Ma are assigned 
to 


Cee es Xt actge lor XU 


Ly lq 


respectively. Here 2 xij is the sum of the xij’s corresponding to the dij’s_ belonging 
L 


k q 
to Lx. The priority factors Mx’s are all positive and are such that the expression = 
=0 


ak Mx has the same sign as the nonzero a« with the smallest subscript present in it 
irrespective of the values of other Li’s. Having done all this, the Ist two-objective 
problem of the given two-objective location problem is reduced to the equivalent Ist 
single-objective transportation-type problem of the given two-objective location problem 
seeking to determine xij ~ OG =a mn tom lye + 1) which minimize 


MymM4i, 


Nal ™m™ 
Z1 = Mo[kmyi + (2 bs — Zat)emye + xr SZ sz x15] 
j=1 i=l i=1 j=l 


Mt X xij eich) 


4 


! Me 


k 


subject to the equations (9), (5), (11), (12). This Ist single-objective transportation- 
type problem is amenable to solution by the standard transportation method discussed 
by Hadley!. The minimum value of the total cost Cy comprising the capital ex- 
penditure to be incurred in constructing a new depot at potential site | plus the present 
value of the expenditure to be incurred in total dead mileage over the planning horizon 
is obtained by the coefficient of Mo in the expression on the right-hand side of eqn. (13) 
after it has been minimized. And if Mu is the priority factor with the smallest sub- 
script among all the priority factors assignel to the nonzero sums = xij in the ex- 
pression on the right-hand side of eqn. (13) after it has been minimized, the the minimum 
of the maximum distance D1 among the distances traversed by individual buses from 
the depots to the starting points of their respective routes is givee by the value of the 


dij’s belonging to Lu. 


4. NuMBRICAL EXAMPLE 


Now the above procedure is applied to obtain the optimal solution of a numerical 
problem which is obtained by taking bey1 = 5, C241 2000, cox2 = 2100, koy1 
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=200000, ko,2 = 250000, m = 2,n = 6.N = 20,p = 6,r = 10, s = 2, and assign- 
ing numerical values to all other quantities in the location problem formulated in 
Section 2. The tableau representation of the numerical problem is shown in Table I. 


In this Table, rows with the headings ‘Depot 1’ and ‘Depot 2’ refer to existing 
depots while entries in these rows in the column with the heading of ‘Capacities of 
depots’ refer to their capacities. And rows with the headings ‘Depot (2 + 1)’ and 
‘Depot (2 + 2)’ refer to depots at potential sites while entries in these rows in the 
column with the heading ‘“‘Capacities of depots” refer to their capacities which are 


TABLE I 


Capacities of depots aud buses required at starting points of routes and distances from depots 
to starting point of routes 
sss 


Starting points of Capacities 

routel route2 route3 route4 routes route6 of depots 
Depot 1 2 l me | 5 4 3 25 
Depot 2 1 2 3 5 2 20 
Depor (2 + 1) 3 fe 1 2 a 4 53+ 5—45=11 
Depot (2 + 2) is 2 I 3 2 3 51+5—45=11 
Buses 
required 12 10 8 6 7 8 
_— > 

TBALE IL 


Tableau providing optimal solution of Ist single-objective transportation- 
problem 


type problem of numerical 


Starting points of Spare Capacity 


capacity of depots 





available 
route | route 2 route 3 route 4 route 5 route 6 

74578.8Mo9 37289.4My 1 8644.7Mo 186447.0My 149157.6Mo 111868.2My 0 

Depot | +M, +My; +Me +MI1 +Mz, +M3 25 
10 8 2 5 

37289,4My 74578 8Mo 111868.2Mg 111868.2M, 186447.0Mo T4578.8My 0 

Depot 2 +M;, +M, + M3 + Mz +M;, +M, 20 
12 0 8 

111868.2My 18644.7My 37289.4Mo . 74578.8My) 111868 2Mo 149157.6My 0 
Depot 
(2+ 1) +Mz +Me +Ms +M, +My; +M, 1] 

6 5 
Bucs 

required 12 10 8 


So 
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identical and equal to 5] + 5— 45 = 11. For easy reference. the term cell shall be 
used to denote the space in a column with the subheading of a route against a depot. 
The entry of all the cells (i, 7) depicts the units of distance dij from depot ito the 
starting point of route j. The two objective functions of the numerical problem which 
are sought to be minimized are 


2 6 2 6 
C= = doyt[koyt + Il co + ZT coyni x(243) + 24 ie Cij Xtj 
I=1 j=1 i=1 j= 


...(14) 
where 
cij = 2 (365) (6) dij = (1 + 10/100)-§ G@ = 1, ....2 + 2; 7 = 1, .... 6) 
ha Ato) 
and 
D = max {dij : x43 > 0 (i = 1, ...,.2 +2; j=1,...,6)}. ELL) 


For the Ist two-objective problem of the numerical problem, we find g = 6. The 
subsets forming partition of the set {dij :i = 1, 2.2 + 1; j = 1, ..., 6} are as follows: 


Ly = {dy4, do5}, Le = {d15, dse}, Ls = {416, d23, doa, 431, 435} 
La = {1, doz. dog, dz4}, L5 = (12, doi, d33}, Le = (dis, 432} 


The di;’s belonging to Z1, Lo, L3, L4, L5, Le have numerical values 5, 4, 3, 2, 
1, .5 respectively. The objective function of the Ist single-objective transportation- 
type problem of the numerica) probiem is given by 


{ Mo (222000 + 74578.8x11 + 37289.4x12 + 18644.7x13 : 
| +186447.0x14 + 149157.6x15 + 111868.2x16 + 37289.4x21 : 
: + 74578.8x92 + 111868.2x23 + 111868.2x24 | 
| + 186447.0x25 + 74578.8x26 + 111868.2x31 + 18644.7x32 | 
: 4+ 37289.4x93 + 74578.8x34 + 111868.2x35 149157.6x36 
: 

| 


Zi = 
"4 


| 
+ My (xia + x25) + Me (x15 + X36) Ms (X16 + X24 | 
| + x31 + x35) + Ma (x11 + x22 + x26 + X34) + Ms (x12 | 
l + xo1 + x33) + Me (x13 + x32) J 

(17) 

Applying the standard transportation method and remembering that the priority 


6 : 
factors Mx’s are such that the expreesion & Lx Mx has the same sign as the nonzero 
k=0 


Lx with the smallest subscript present in it irrespective of the values of other Le Ss; 2 
optimal basic feasible of the Ist single-objective transportation-type problem is obtain- 
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ed and is shown in Table II. Entry in the upper half of the cell in Table II depicts the 
cost associated with the variables. Values of the basic variables of the optimal 
solution are encircled in the tableau. The optimal value of the objective function Z1 
is given by 

Z, = 3093283.8Mo + 2M2z + 5Mz3 + 14Mq + 22Ms5 + 8Me6 .-.(18) 


An optimal basic feasible solution of the 2nd single-objective transportation- 
type problem of the numerical problem is obtained proceeding exactly inthe same 
way as done to obtain the optimal solution of the Ist single-objective transportation- 
type problem of the numerical problem. The optimal value of the objective function 
Zz2 of the 2nd single-objective transportation-type problem is given by 


Zz = 3107094.4Mo + 6Mz3 + 17M5 + 20Me6 + 8M7 sue b9} 


Now among the two single-objective transportation-type problems of the nume- 
rical problem, the objective function of the Ist one has minimum value. So the site 1 is 
selected for construction of a new depot and the optimal solution of the numerical pro- 
blem is provided by the optimal solution of its Ist single-objective transportation-type 
problem. For ready reference, the optimal solution of the numerical problem is 
shown in Table III. 


TABLE Ill 


Optimal solution of numerical problem 


Optimal site Values of basic Minimal value of Minimum of 
for new depot variables of optimal expenditure in maximum of dead 
solution construction plus mileage of individual 
present value of buses 
expenditure in 
total dead 
mileage 
l/=1 x12 = 10, A jg) 8x15 = 2 C = 3093283.8 D= 4 
x17 = J, x2] = 12, x2 = 0. 
Xe = 8, tans x35 = 5 
ll 
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A NOTE ON NORMED NEAR-ALGEBRAS 


T. SRINIVAS AND K. YUGANDHAR 


Department of Mathematics, Kakatiya University, Warangal 506009 (A. P.) 
(Received 29 March 1988) 


In this paper an attempt has been made to extend the theory of Algebras to 
Normed near-algebras. For this purpose we are led to a definition of a 
Normed near-algebra which differs from previous definition. A technique 
has been developed to find the regular elements in the Normed near-algebras 
and proved that the set of all regulas and quasi-regular elements is open. 
Further it has been shown that the concepts of topological divisors of zero 
and the spectrum of an element can be carried to Normed near-algebras 
also. 


§1. Near-algebras were studied by Yamamuro®, Brown? and others (cf. Pilz’). 
In order to extend the theory of Normed algebras to Normed near-algebras we have 
defined the Normed near-algebra with an additional condition. 


Definition 1—A (right), near-algebra B over a field F is a linear space over F on 
which a multiplication is defined such that(i) B formsa semigroup under multiplica- 
tion, (ii) multiplication is right distributive with respect to addition, (iii) «(xy) =(«x)y 
for all x, y € Banda € F. 


Definition 2— A near-algebra B over the real or complex numbers is called a 
Normed near-algebra provided that there is associated with each x € Ba real num- 
ber || x ||, called the norm of x, with the following properties : 


(i) || x || > O and || x || = 0 if and only if x = 0 (the additive identity of B), 
(ii) Ix + yl Sxl +1 yl, for all x, y € B, 
(iii) || «x || — | « | || x |], «a real or complex number, 
(iv) |! xyll <x I I» |, for allx,y € B, 
(v) ||xy — xz <x ly — z\, forallx,y,z€ 8B, 
(vi) If B has an identity e, then || e || = 1. 


Through out this paper the property (v) plays a crucial role in proving most of 
the theorems. It follows from the definition (2) that all the operations are continuous. 


Through out this paper we assume that 4 is a normed near-algebra such that 
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(i) Bis complete with respect to the norm defined as above, 
(ii) B has an identity e, 
(iii) Forevery x € B, 0.x = x.0 = 0. 


Example— Let S be a Banach space over the field F of real or complex num- 
bers. Then the set B (S, S) of all Lipschitz continuous functions x : S — S such that 
x (0) = 0, and having the Lipschitz norm and with the composite map x (y (s)) 
= (xy) (s) of S onto S, for allx, y € B(S,S) ands € S, constitute a normed near- 
algebra‘, 


Definition 3— An element r € B is said to be left (right) regular if there exists 
an element s € B such that sr = e (rs = e). 


The element s is called left (right) inverse forr. An element which is both left 
and right regular is called regular element. An elemert which is not (left, right) 
regular is called (left, right) singular. 


Through out this we denote the set of all regular elements by G and the set of 
all singular elements by S. 


The following Lemma is very useful in proving the main theorem. 

Lemma |— Let B be normed near-algebra, then every element r € B such that 
|e —r || < 1 is regular. 

Proor : Let x = e — rand yo = e. Define the sequence {yn} “~, inductively by 


yng1 = Yo Xyvn. 
Then 


l| Yn+1 — Yn|| = || xyn — xyn-1 || 
S || x Il ll yn — yn-1 |) 
< || x i" lyr — yo ll 
== || x *"2 
and hence || yn+p — yn || < || x |"? + ... + || x |l"*1 for any integer p. 


Therefore {jn} is a Cauchy’s sequence in B since || x || < 1. As Bis complete, 
there exists a »y € Bsuch that Lim yx = y. From the relation yn.i = Yo + xyn, we 


n->oo 


gety = e+ xy =e + (e — r) y, that is ry = e, 
Next consider 
ll ynr — e|| = || xyn-t.r— x] 


<x Il ll year —el| 


Fe 
<= lx er. 
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Therefore, || yr — e || = Lim || ynr — e || = 0, that is yr = e. 

n>o 
Hence r is regular. 


Remark 1: Here the inverse of r is given by r-! = Lim yn, and it can be 
n—->oo 


proved by using the same construction for yn that 


1 


r-1|| << ——_____. 


Lemma 2— Let G1 denote the set of all left regular elements of B. If x € Gi 


then any element z € B satisfying || x — z || < || y ||"! (where y is a left inverse of x) 
belongs to Gi. 


Proor : || e — yz || = ll yx — yz || 
<|lyililx — zl 
a 


This implies by the Lemma (1) that yz © G and hence z € Gi, which proves the 
Lemma. 


Similarly we can state : 


Lemma 3— Let Gi denote the set of all right regular elements of B. If x € Gi 
such that the right inverse of x is y then any element z € B satisfying ||x—z|| < ||yl|-1, 
belongs to Gr. 


Now we will prove the main theorem of this paper. 
Theorem 1— The set of all regular elements G is an open set in B. 
Proor: Letx € G = Gif Gr. 
Consider the open ball 
S (x) = {z € Bill z — x || < 1/l] x1 Ib. 


Then for any z € S (x), by Lemma (2) and Lemma (3), z is regular. Hence 
S (x) C G. Therefore G is open. 


Remark 2: Infact it can be proved that Gz and Gr are open and hence 
G = Gi (1) Gr is open. 


Theorem 2— Let B be the normed near-algebra. The mapping r > r-lisa 
homeomorphism of G onto G. 


Proor : It is sufficient to show that the mapping r > r~! is continuous. 
Let x, r € G such that 
Wx — rj] < 1/2 |r? I. 
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Then 
he — xr-2 |] = [)rr-t — xr-* | 
iin — x ie 
< lf2, 
Hence xr-! € G and || (xr~4)-1 || < 2, by Remark 1. 
Consequently, 
I] x-2 |] = [rod rx | 
< |] rll rx-7 |] 
= |r Ulil@rty?| 
<2 || r77 Il. 
Therefore 
phic Need seria eo eer ae 
<|}x2ix—r ie tl 
<2] r2|P Il x — ri. 


This shows that the mapping r—>r-! is continuous. Thus we have r —~r-! is 
homeomorphism. 


Remark 3: The above theorem implies that the set of all regular elements G is 
a topological group. 


The following theorem can be proved by using the same techniques as in 
algebras. 


Theorem 3— Let B be a normed near-algebra and {rn} a sequence of left (right) 
regular elements of B which converges to anelementr € B. If sn is a left (right) 
inverse for rn and if {sn} is a bounded sequence, then r is also a left (right) regular, 


Proor : Let {rn} be a sequence such thatrn > rin B. Let {sn} be a bounded 
sequence that is sup || sx || = M < co and snis a left inverse of rn. 
n 


As rn > r we get, for « = 1/M > 0 there exists a number 79 > 1 such that 


|rn —r|| < 1(M + 1) forn > no. 
Now consider 


|e — snyr || = || Sng Ty — Sno r || 
< [I sng lll ray — r | 


(equation continued on p. 437) 
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< || sn, "(M+ 1) 


Gedc 


This implies that sn) r € G, thatisr € Gi. 


§2, Definition 4— Let B bea normed near-algebra. An elementr € B is said 
to be quasi-regular if e — r is left invertible. (cf. Beidlemen and Cox!). 


Lemma 4— Let B bea normed near-algebra, then every element x € B such 
that || x || < | is quasi-regular. 


Proor: This is obvious, since || x || < 1 implies that || e — (e — x) Ms 
Hence by Lemma 1, e — x is invertible, that is x is quasi-regular. 


Theorem 4— Let Q denote the set of all quasi-regular elements of B. Then Q is 
open in B. 


Proor: We know that Gz is anopen set in B. But an element r € Q<+(e—1r) 
€ Gi. Hence Q is open in B. 


Definition S— An element z in the normed near-algebra B is called a left (right) 
topological divisor of zero provided there exists a sequence {zn} in B such that 
| zm || = 1, for all n and zzn > 0 (zn z > 0). 


We call the element which is either a left or right topological divisor of zero 
as a topological divisor of zero. We denote this set by Z and the set of all left (right) 
topological divisors of zero in B by Zz (Zr). 


We now state the following theorems without proofs since they are very easy 
and can be proved using the same techniques as in algebras. 


Theorem 5— The set of all topological divisors of zero is a subset of the set of 
all singular elements. In other words in the normed near-algebra B every topological 
divisor of zero is singular. 


Theorem 6— The boundary of S in B is a subset of Z. 


Theorem 7—(a) Every element z for which (e —z)is a topological divisor of 
zero is quasi-singular, that is not quasi-regular. 


(b) Gr C Hi, Gr C Hy and G C H, where Hi, Hr and H are the sets that are 
complements of the sets G7, Gr and G respectively. 


(c) G= Gi Ar = Gr M1. 
(d) Ge i8e GC Zr. 


Definition 6—Let Ba complex normed near-algebra and let x be any element of 
B. Let 50 (x), be the set of all non-zero complex numbers such that e — xA~le is 
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singular. Put o (x) = so (x), if x is not singular and o (x) = a9 (x) U {0}, if x is 
singular. 

Then oc (x) is defined as the spectrum of x. 

The complement of « (x) is denoted by P (x), 

Theorem 8— sc (x) is bounded and closed in B. 


Proor: Let A € 59 (x). By the Lemma 1, we get || e — (e — xA-le)|| ¢ 1. 
That is || xA-1 e || ¢ 1, which implies that || x || > | A | . Hence o (x) is bounded. 


The mapping A > xA-le is continuous mapping from C — {0} to B. Hence 
Theorem |, gives that og (x) is closed and hence ¢ (x) is closed. Thus the proof is 
complete. 


But it is an open problem to settle whether o (x) is non-empty or not when x is 
not singular. 
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Here we prove the existence of the classical near-ring of right quotients of 
the sub near-ring D of the distributive elements of a prime Goldie Abelian 
near-ring K which is such that some power of each element of K belongs to 
D and further its non nilpotent elements are distributive. 


§1. We have already introduced the notion of a Goldie near-ring and have 
established some results on prime and semiprime Goldie near-rings, proved Goldie 
Theorem analogue for Abelian Goldie near-rings in which non nilpotent elements are 
distributive!. 


In this paper we discuss a radical Goldie near-ring. If A is a subnear-ring of a 
near-ring K then K is radical over A (or K is A-radical) if a power of each element in 
K lies in A. In an Abelian near-ring in which non nilpotent elements are distributive? 
K is D-radical, where D is the subnear-ring of distributive elements of K. Here we 
establish that if K isa prime Goldie Abelian near-ring which is radical over the sub- 
near-ring D of distributive elements of K and further its non nilpotent elements are 
distributive, then D has a classical near-ring of right quotients. Unless otherwise 
specified a near-ring K will contain unity and a.0 = 0 for alla € K. 


§2. A countable ordered family {Ai, A, ...,} of subsets of a near-ring K is an 
independent family if for all n EN, A O( 3 Ar) = 0, were 1 (i < nand 
kFi 


beat Gens 


A right (left) K-subset 4 of (right) near-ring K is a subset of K such that 
ak (ka) € Afora€ A,k E K. If A is a right and a left K-subset of Kthen Aisa 
K-subset of K. A is a right (left) essential K-subset of K if for any non zero right (left) 
K-subset of K has a nonzero intersection with A. A near-ring K is semiprime if K has 
no nonzero nilpotent K-subset, K is prime if O is a prime K-subset (i. e. for any two 
K-subsets A, Bof K, AB =0 > 4 = 0 or B = 0). 


A Goldie neart-ring is a (right) near-ring K with a.c.c. on right annihilator K- 
subsets and which has no infinite independent family of right K-subsets of it. 


In an Abelian near-ring K, the subset D of distributive elements of K is a sub- 
near-ring of K (it is actually a ring). The following two Lemmas are easy. 
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Lemma 2.1.1— Let K be an Abelian near ring and S bea set of distributive 
elements of K. Then the right annihilator K-subset rs (S) is a right ideal of K. 


Lemma 2.1.2— If A is a subnear-ring of a near-ring K and S C A, then 
ra(S) = re (S) 1 A. 


Lemma 2.2.1— Let K be a near-ring. Suppose k € K is such an element that 
for each x € K there exists n € Z* (depending on x) satisfying the condition 
kxn k = 0. (In Lemma 2.2.7 we see that such a condition is justified). Then there 
exists a nonzero element a of K with a2 = 0 and satisfying the same hypothesis as k. 


Proor : Assume x = k. Then k* = 0 for some « > 2. For k?2 = 0 we get the 
required result on taking k = a. And if k2 0 let « be minimal such that k* = 0. 
Then a = k*-1 (+ 0) satisfies the required hypothesis. 


Lemma 2.2.2— Let K be a semiprime Goldie near-ring and k € K, k # 0 be as 
in Lemma 2.2.1. 


Then there exists a € K with a2 = 0 such that a satisfies same hypothesis as k 
and for any m € K with m= = 0 we have ama = 0. 


Proor; By above Lemma there exists a € Ka 0, a* = O satisfying the same 
hypothesis as k. Let m € K be such that m2? =0. Then by given hypothesis there 
exists « > 1 such that for any x € K we have 


a(m + maxa)* a = 0, or a (m + maxa)*"1(m + maxa)a = 0 


or, a (m + maxa)*~! (ma) = 0 (since right distributivity holds in K and a? = 0). 
Using the fact that m? = 0 and repeating in like way we get finally (xama)* = 0. 


Thus, Kama is a nil left K-subset of K. And K being semi-prime we therefore get 
that ama = O(I1E &K),. 


Lemma 2.2.3— Let K be as in Lemma 2.2.2. Then there exists a € K,a + 0, 
a* = 0 such that a satisfies the same hypothesis as k in the above Lemma and for any 
P,gd © K, pq = O implies paq = 0. 


Proor: From above Lemma we get a € K such that a2 = 0,a ~ Oanda 
satisfies same hypothesis as k and for any m € K with m2 = 0 we have ama = 0. 


Now let p,g = Kand pg = 0. Then for any x € K, (gxp)2 = 0 which gives 
a (9xP) a = 0 for any x € K. Thus (Kpaq)2 = 0 from which Kpaq = 0 follows for 
K is semiprime Goldie. Therefore pag = 0 (since 1 € K). 


Lemma 2.2.4— Let K be a semiprime Goldie near-ring and k € K be such that 


— each x € K, there exists ann € Zt (depending on x) such that kx™ k = 0 Then 
== 0, ' 


Proor : If k ~ 0, then by Lemma 2.2.3 we geta Ee K,a540, a2'= 


a satisfies the hypothesis as in the Lemma. O such that 


Now for any x € K we getn € Z+ such 
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thatax"a = 0. Ifm = 1 then axa = 0. Therefore (xa)? = 0. And ifn > 1 then 
ax" a = 0 gives (ax"-1) a(xa) = 0. By the hypothesis satisfied by a we get (ax”~1) 
a (xa) = 0. And this gives ax”~2 (xa)? = 0. In like manner we get ax"-3 (xa)? = 0. 
Finally we get (xa)"*! = 0. Thusin any case we get (xa)"*1=0. So Ka is a 
nilpotent left K-subset of K. And K being semiprime it therefore follows that Ka = 0 
which gives a = 0, a contradiction. Hence k = 0. 


Lemma 2.2.5— Let K bea prime Goldie near-ring. Suppose a,b € K,a #0 
and for each x € K there exists ann € Z* (depending on x) such that ax" b = 0. 
Then 5b = 0. 


Proor: Letp = {fy € Kl y x»@) b = 0, x € K}. Since for any k € K, 
kyxn(*) b = 0 is a left K-subset of K. And from Lemma 2.2.4 (by) x"@ (by) = 0 
gives by = 0. Thus b P = 0 which gives (KbK) (PK) € Kb? K(= 0). And since 
pK + 0 and K is prime it therefore follows that KbK = 0. And hence b = 0. 


Lemma 2.2.6— Let K bea prime Goldie near-ring and A be any subnear-ring 
K. If Kis A-radical, then 


(i) for any A-subset J, J of A with JJ = 0, J] = 0; we get J = 0; 
(ii) if a € A andra (a) = 0 then r« (a) = 0. 


Proor: (i) Leta € J,a~O0andb € J; K being A-radical for x € K there 
exists © Z+ such that x” € A. Now ax" b € JAJ C IJ (= 0). Thus ax” b = 0, 
a + 0. Therefore by Lemma 2.2.5 we get b = 0, So J = 0. 


(i) If rx (a) #0, let x € re (a), x #0. Then for some n € Z* we have 
xm € A and ax” = 0. So x” € ra (a) which means that x is nilpotent. Thus r« (a) 
is a nil right K-subset of K. And K being prime it therefore follows that r« (a) = 0, 


Lemma 2.2.7— Let K be a semiprime Goldie near-ring and A a subnear-ring 
of K such that K is A-radical. Then 


(i) A cannot have any nonzero nil right (left) A-subset; 
(ii) if for ay, ag € A we have aj A (\ a2 A = 0, then a1 K M ag K = 0. 


Proor: First we show that A cannot have nilpotent right A-subset. Let P be a 
nilpotent right A-subset of A. Without loss of generality we may assume P? = 0 and 
leta € P,x © K. We getn € Z* such that x" € A. And then ax®a€ pAPC 
pp (= 0). Thus for any x € K there exists m © Z* such that ax" a = 0. And by 
Lemma 2.2.4 we get a = 0. Therefore Pp = 0. Now we show that A can not have 
nonzero nil right A subset. For this we show that a satisfies the a.c.c. for right 
annihilators. 


Consider an ascending chain 


ra(Si) Cra (So) C ... of right annihilator; where Rit San s-.5'G A. Then 
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re (la (ra (S1))) © rx (la (ra (S2))) C .... Since K is Goldie, we therefore get t © Zt 
such that 

rx (la (ra (St))) = r« (14 (ra (St41))) = ... 
or, 

AN (rg (La (ra (St))) = A 1 (re (Ia (ra (St41))) = «.-. 


And because of Lemma 2.1.2 we get 


ra (la (ra (St))) = 14 (la (ra (Sti) = «-- 
Hence 
ra (St) = ra (St,1) = «-- - 

(ii) Consider the right A-subset P = a; K 1) a2 A. Letx € P. Then x = ak, 
= ay ko where ki € K, ko € A. Since K is A-radical there exists > 1 such that 
(ky ay)" € A. Since a1, ko € A we get a (ky ay)” = (ay ky)" ay = (ag ko)" a € 
a, A (\ a2 A (= 0). 

Therefore 
ay (ky a1)" = 0. 
Hence 
xmt1 = (ay ky) (a1 ki)™ = ay (ky a1)" ky = 0. 


Thus P is a nil right A-subset of A. Therefore P = 0. Soa; KQ a2 A = O. 


Now let x € a1 k (1 a2 K and suppose x = aj ty = a2 te, where ty, tg € K. 
Since K is A-radical we get m € Z* such that (tog a2)" € A, Then (ay t1)™ ae 
= (a2 t2)™ az = ag (te a2)™ E ay K () ag A(= 0). Soxm+l = (ay t)™+1 = O, Thus 


a, K ()\ a2 K is a nil right K-subset of K. K being semiprime, we therefore get 
a,K (\aeK= 0. 


The following Lemmas have been proved in Chowdhury!, 


Lemma 2.2.8 : For any right essential K-subset “A of a near-ring K, ifa € A 
then 


a1lA={x€ K|ax€ A} 
is a right essential K-subset of K (Corollary 2.3.3 of Chowdhury}), 


Lemma 2.2.9— Let K be a Goldie near-ring whose non-nilpotent elements are 
distributive. If x € K be such that r (x) = 0, then xK is a right essential K- 


subset of 
K (Lemma 2.4.1 of Chowdhury!), 


Lemma 2.2.10— A semiprime right Goldie near-ring is right non-singular 
(Lemma 2.4.3 of Chowdkury!), 


Lemma 2.2.11— If Kis a semi prime Goldie near-ring where non nilpotent 
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elements are distributive then every right essential K-subset of K contains a regular 


element. (Lemma 2.4.4 of Chowdhury!). 


Lemma 2.2.12— Let K be a semiprime Abelian Goldie near-ing such that its 
non-nilpotent elements are distributive. Then, K-satisfies the d.c.c. for right annihilator 
ideals (Lemma 2.4.5 of Chowdhury!). 


Lemma 2.3.1—Let K be a prime Goldie Abelian near-ring which is such that its 
non-nilpotent elements are distributive. If Dis the sub near ring of distributive ele- 
ments of K, then D satisfies the D.C.D. on right annihilators. 


Proor : If S C D, then rp (S) = r« (S) (1 D. 


Now let rp (S1) D rp (S2) D ..., where Sj, Sa, ... C Dbea decending chain of 
right annihilators in D, Then 
r« (Ip (rp (S1))) 2 rx (Ip (rp (S2))) 2 --- 


Here rx (Ip (rp (S1))), ... etc. are right ideals. Therefore by Lemma 2.2.12, these 
exists t © Z* such that 


rx (Ip (rp (St))) = rx (Ip (rp (St+1))) = --- 
Therefore D ( rx (Ip (rp (St) = DO re (I (rp (St41))) = --. 
And by Lemme 2.1.2 we get rp (1p (rp (S:))) = rp (Ip (rp (St,1))) which gives 
rp (St) = rp (Sty) = --- 
Thus D satisfies the D.C.C. on right annihilators. 


Lemma 2.3.2— Let K be as above. If for some a € D, we have rp (a) # O then 
aD is not a right essential D-subset of D. 


Proor: If aK <e K, then by Lemma 2.2.11 re (a) = 0 and therefore rp (a) 
= rg (a) (\\ D = 0, acontradiction. So aK Xe K. If possible let aD <e D and con- 


sider any nonzero right K-subset C of K. Then C is not nil and let c bea non-nil- 
potent element of C. By hypothesis c € D. Since 1€D,cD#0.SoaDN cDF 0. 
Now 
aK 7. CDaKNcKDaDNcDF 0. 
Thus aK <e K,acontradiction. Hence aDe D. 


Lemma 2.3.3— Let K and D be as above and Kis D-radical. Then a right D- 
subgroup A of D which is also a right essential D-subset of D contain a regular 


element. 
Proor : Choose a € A so that rp (a) is minimal in D (Lemma 2.3.1). We 


claim that rp (a) = 0. 
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If rp (a) ~ 0 then by Lemma 2.3.2 aD te D. So there is a nonzero right D- 
subset J of D such that aD () J = 0. Since A <e D, 


AN JF 0. Write J’ = AM J. 
Then 

aDNJ=aDC(ANJ)=(@DNINA=O. 
If x € J’ andd € rp (x + a), then xd + ad = 0 
Or, 

ad = — xd = y(—d) (since x € D) 

= 0 (for aD (\ J = 0). 

Thus d € rp (a) (\ rp (x). Therefore rp (x + a) € rp (a) ( rp (x) C rp (a). 
And x + a € A (since A is a subgroup of D). 


Hence we get, by minimality of rp (a), rv (x + a) = rp (a). 


So rp (x + a) = rp (a) (\ rp (x) = rp (a) which gives rp (a) C rp (x). Hence 
x rp (a) = 0. Therefor we get J’ rp (a) = 0, Thus (DJ’) (Drp (a)) € DJ’ rp (a) (=0). 
So Drp (a) = 0 for DJ’ # 0 (since 1 ED) (Lemma 2.2.6 (i)). It follows that rp (a)=0. 
Therefore rc (2) = 0 (Lemma 2.2.6 (ii)). New letc € lx (a). Then ca = 0 which 
gives caK = 0. By Lemma 2.2.9 aK <e Kand therefore c € Z (K). And K being 


right non singular (Lemma 2.2.10), we get c = 0. Thus lx (a) = 0. It follows that a 
is regular in D. 


§3. We now prove the main result. 


Theorem 3.1: Let K be a prime Goldie Abelian near-ring which is D-radical 


and whose non nilpotent elements are distributive. Then D has a classical near-ring 
of right quotients. 


Proor: Let a,b € D and a be regular in D, Thenrp (a) = 0. By Lemma 
2.2.6 (ii) we get rx (a) = 0. So aK €e K (Lemma 2.2.9). Now if aD <e D, then there 
exists a non zero right D-subset P of D such that aD A P= 0. Letx € P, x0, 
we then get aD (} xD C aD (\ P (= 0) oraD Q xD = 0 which gives aK () xK = 0 


(Lemma 2.2.7 (ii)) put aK <e K then gives xK = 0 whichis a contradiction (since 
1€ K). SoaD <e D. 


Now let A = {x € D| bx € aD}. Then by Lemma 2.2.8, AX <e D. And dA is 
a subgroup of D (Since b is distributive and aD isa right ideal of D). Thus A is a 


right D-subgroup of D. So by Lemma 2.3.3, A contains are 


gular element, Say a’, 
Hence ba’ = 


ab’, for some b' € D. Thus the right Ore condition is satisfied with 


respect to the set of regular elements in D. Therefore by Theorem 3.31 D has a 
classical near ring of right quotients. 
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The purpose of the present work is to derive a few results which improve 
upon some earlier finding on stability analysis of solutions of ordinary differ- 
ential systems. Firstly, wetreata most general differential system with 
impulsive perturbations and study the stability of the solutions of such a 
system. The results thus derived unify and improve the results of Raghvendra 
and Rao‘ and Strauss and Yorke simultaneously. Further we also treat the 
Stability problem of a weakly linear differential system which is new and infact 
is an extension of result of Raghavendra and Rao. One final result is derived 
on necessary and sufficient condition tor stability of specific nonlinear differen- 
tial system. 


1. INTRODUCTION 


Raghavendra and Rao? investigated some Stability properties of solutions of 
ordinary differential systems with respect to impulsive perturbation. In fact slight 
modification of Gronwall-Bellman type integral inequalities and their extentions, which 
has been given by them, is the main tool in considering the stability theorem for solu- 
tions of ordinary differential equations containing measures. 


In the present work we try to deal with most general type of differential system 
analogous to Strauss and Yorke® with a difference that we now introduce impulsive 
perturbation in the system considered by him. Subsequently the Stability theorem thus 
developed extends in theory as purported by Raghavendra and Rao‘, The next section 
is devoted to preliminaries and basic lemmas quoted from Raghavendra and Rao’ and 


Strauss and Yorke5 which we have used later on in the sequel. Subsequent sections 
deal with the main results of the present work. 


2. PRELIMINARIES AND Basic LEMMAS 


: n 
For any vector x € R®, the euclidean space of dimension n,let|x| = 3 | x¢ 
i=] 


By C [E, R”] we denote the Class of continuous mappings from E into Rr. 
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Let J = [0, co]. Most general type of differential system with impulsive pertur- 
bation is given as follows. 


Dx = f(t, x) + g(t, x) Dv + h(t) note 1) 


where x € R”, Dy denotes the distributional derivative of the function v,f,g € C 
[J x R", RX} andv:J— Ris a function of bounded variation and it is right 
continuous on J. Dv can also be identified with a Stieltjes measure and in fact 
possesses the effect of the sudden change of the state of the system at the points of 
discontinuity of u. 


Equations (2.1) may be regarded as a perturbed system of ordinary differential 
equation 
Dx = f(t, x) ...(2.2) 


where the perturbation g (t, x) Dv is impulsive and h(t) satisfies certain smoothness 
condition of the type as discussed in Strauss and Yorke’, to be stated later on in the 
sequel. We also assume that g (f, x), f(t, x) and A(t) satisfy certain smoothness 
conditions sufficient to guarantee the existence and the uniqueness of the solutions of 


(2.1). 
We assume the following conditions : 
(Gy) There exists « > 0 such that if |x| <a, then | g(t, x) | < y (@) forall 2 0, 


where 
G (t) = te y (s) dv (s) -> 0 as t + ©© ARR) 
t 


where y € C[J, R,]. 
(Gz) There exists a continuous, non-increasing function H (ft) satisfying : 


Lim H(t)=0 


t-o 


t 

such that | f h(s) ds | < H (to) for every 
t 
0 


0<to St<stotl. 


We shall prove that if g (t, x) satisfies (G1) and / (t) satisfies (G2) then there exists 
Ty) > 0 and 59 > O such that if fo 2 To and | xo | < 8, the solution F (1, to, xo) of 
(2.1) approaches zero as ! > ©. In particular if x = 0 is a solution of (2.1) then it is 
uniformly asymptotically stable. In fact on the other hand if g (t, x) = 0 and A (t) 
does not satisfy (G2), then no solution of (2.2) tends to zero as f > %. 


In the case f(t, x) = Ax, where all the characteristic roots of A have negative 
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real parts, we derive some stability results for the so called weakly nonlinear differen- 
tial system of the form 
Dx = Ax + F(t, x) + g(t, x) Dv + A(t) .-(2.4) 


we now give some definition in the next paragraph. 
Definition 2.1—The null solution of (2.1) is said to be uniformly asymptotically 
stable if the following two conditions hold : 
(i) forevery « > 0, there exists a 5 = 5 (ec) > 0 and 
+ = +(e) > O such that | y(t, to, x0) | < «, 
1 > to > + (€) provided that | xo | < 3. 
(ii) for every » > 0, there exists positive numbers 
8o, to and T = T(n) such that | y (t, to, xo) | <7 
t>to + T, to > to provided | xo | < ‘to. 


In the next section we would extend the proof of the Theorem 3.1 of Raghavendra 
and Rao‘ to get our result and the proof thus given does not make use of the 
Lyapunov function as also done in the case of Theorem 5.1 of Strauss and Yorke®. 
Before we prove our result we give some lemmas which are basically the extensions of 
the lemmas which appear in Strauss and Yorke® and we also state the Gronwall type 
of inequalities quoted from Strauss and Yorke® which has been used further in our 
result. We also remark here that the generalized inequality of Raghavendra and Rao4 
is not needed in the proof of our result, instead the more traditional inequality such as 
the one given in Lemma 2.4 below suffices for the derivation of Theorem 3.1. 


t t 
Lemma 2.1— J Y (s) dv(s) < J G (s) ds for all t > to > 1. 
0 0 


t t 
Lemma 2,2— J evs y (s) dv(s) < J : ev(stl) G(s) ds foallo > 0 ¢ > to > 1. 
0 o- 


t 
Lemma 2.3—lim  e-? ! ers y (s) dv(s) = 0 forallo > 0. 


too 


Lemma 2.4—Let » (t) and p(t) be non-negative and continuous for ¢ > fo, let 
C > 0, k = 0, and let 


t 
y@)<C+ J [k y (s) + p(s)] ds 
then 


zm t 

y(t) << CeO) 4 FD (5) ek (t-2) ds. 
t 
0 
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The following theorem gives the sufficient conditions for uniform asymptotic stability 
of (2.1) with respect to (2.2). Let f(t, 0) = 0 forall t > 0. 


Theorem 3.1—Let the null solution of (2.2) be uniformly asymptotically stable. 
Assume that : 


(i) f satisfies the uniform Lipschitz continuous condition of the following type. 
if) —Jiby a Lix—yl, jel, tyl <a 
where a is a positive real number. 
(ii) g(t, x) satisfies the condition as in (G1) and (2.3) holds. 


Then there exists To > 0 and 59 > 0 such that fo > To and | xo | < 40, the solution 
y (t, to, Xo) of (2.1) satisfies | y (t, fo, xo) | > 0 ast > oe. In particular ifg (t,0) = 0 
then the null solution of (2.1) is eventually asyptotically stable. 


Proor : The solutions and constants corresponding to system (2.2) shall be 
starred, those for (2.1) shall not. Now let 


QO (t) = Sup {G(s): t-—1E5< oo}. 
Then it is easy to see that 
Q(t) 4,0 as.t > ©o. 


Invoking lemma 2.1 we get 
t t 
f y(s)dv(s) < J G(s) ds < Q (to) (t — to + 1) 
to Ae! 


ift > to > 1. Also we have from condition (G2) 


ty+1 


t 
if h(s)ds| <| fi h(s) ds \ + ... --. + | { &A(s)ds| 
to to to+™ 

mH Afb) oF + H(to + m) 


< H (to) (t — to + 1). 
Put B(t) = Q(t) + H(t). Letto > ! and | xo | Gv. 
If y (t, to, Xo) is a solution of (2.1) and if 
| y(t, to, Xo) | < yon (to, to + t] for some y > 0, then 


| y(t, to, x0) — Y* (ts 0, Xo) | 
t 
— | x0 os j f(s; y(s, tO, xo)) ds 37 J g (s, vy (s, 10, xo)) dy (s) 
"0 


0 
(equation continued on p. 450) 
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t 
“t- i h (s) ds — xo — J f (s, ¥* (s, to, Xo)) ds | 
t 0 


t 
<J Liy@—yG) lds + Sy @d +1 J A(ds| 
to 0 0 
<t LiyG)=y* @) | art BU) eo 
t 
0 


From Lemma 2.4, 


| y (t, to, xo) rat vs (t, to, Xo) | 
t 
€& Bi(to) eX) + 5 B (to) e4(¢-8) ds 
t 
0 


S B (to) e& + B (to) ef (t — to) 
= eM (1 +t) B (to). 


We may assume without loss of generality that y < 35 . Next we proceed to show by 
the estimates as derived earlier4’5 that | » (t, to, xo) | < € on every interval [to + mt, 
to + (m + 1) t] and hence on [to, cc). Hence if g(t,0) = Oandh (t) = 0, then it 
is clear that x = 0 is uniformly asymptotically stable. Rest of the proof is completed 


following Raghavendra and Rao4. In fact, given O< y< yy and choosing 
8 (n) = 8* (0/2), O< 38 < 4, 


t (4) = T* (5/2) and Ty (y) so that 

B (Ty) < 8[e# (1 + +) Z}-1 
we finally show that | » (t, to, xo) | < n, forall t > to + 7, 
where T= T(n) = + + Ty. 


4. WEAKLY NON-LINEAR DIFFERENTIAL SysTEM 

We consider a weakly linear differential s 

the present section. The results on stabilit 
4.1 of Raghavendra and Rao4. W 


ystem with impulsive perturbation in 
y derived subsequently generalizes Theorem 
€ consider the following ordinary differential system. 


Dx = Ax + F(t,x) + g(t, x) Dv + h(t) .--(4.1) 
where A isn X n constant matrix and F € C{J x R®, R), The solution ¥(t) of 
(4.1) satisfying y (to) = xo, 19 E J, is given by 


YN =O ~ x0 + [ O— 9) Fy (as 
0 


(equation continued on p. 451) 
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t 
n J © (t — s) g (s, y (s)) dv(s) 
0 


t 
- ! ®(t — s) h(s) ds ...(4.2) 
0 


where © (ft) is the fundamental matrix of the equation x’ = Ax satisfying ® (0) = / 
and ® € C@ (J) where J denotes the unit matrix‘. 


Theorem 4.1—Suppose that 
(i) all characteristic roots of A have negative real parts 
(ii) given « > O, there exists 5 (e), T (<) > 0 such that 
| F(t, x) | <«|x| provided |x| <3(¢) and t > T) 
(iii) The condition (G1) and (Ge) of section 2 hold. 


Then, there exists To > 0 and 59 > O such that for every to > To and | xo | < 4o, 
any solution y(t) = y (ft, fo, xo) of (4.1) satisfies 


lim y(t) = 0. 


a} 


If, in particular (4.1) possesses null solution then the null solution is uniformly 
asymptotically stable. 


Proor : We give a brief outline of the proof. Also note that the technique is 
same as in Raghavendra and Rao4. Because all characteristic roots of A have negative 
real parts, there exists positive constants 4 and ci, such that 


| O(t) | = | e4¢ | < ce for t > O. ...(4.3) 


Assume that 
3 = 3 ((al2) c*) < ¥ 


choose tg and Xo, So that to > T= T ((a/2) C1-4) and | xo | < & <3. So far as 
| y(t)| =|» (t, to, xo) | < 8, for > fo we have 


t 
ly (t)| <c, e*'o | xo | + j (a/2) e-4¢¢ 8) | y(s) | ds 
0 


t t 
+ f cy e-a(t-8) y (s) dv(s) + | J e~4(t-8) h(s) ds | 
vi) 0 
which implies that 


t 
Ly (t) | ett Sere | x01 + a/2 f em | x) | ds 
0 
(equation continued on p. 452) 
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t t 
+ f cy e% x (s) dv(s) + J e%h (s) | ds. ...(4.4) 
i t 


0 0 
Applying Lemma 2.4, we have from (4.4) 


t 
| y (t) ett < ¢, eM | x9 | ef) 4 J ex ea eal2('-#) y (s) dv (s) 
0 


t 

4+ f eti2(t-8) eas | h(s) | ds. 
t 

0 


From which it follows that 


t 
Ly(t) | <1 e *?0) | xo | + 1 f e 22-8) ~ (s) dv (s) 
t 
0 


t 

+ f e7412(-s) | h(s) | ds. 
t 
0 


By choosing 59 and To suitably and using the properties of the function Q (rt) and H (1) 
we can proceed as Raghavendra and Rao‘ to show that for t > To and | xo | < 3, 
we have | y(t) | < 5, which implies the existence of y(t, to, xo) on the interval 


[to, ce). For t > To > 1, the fact that 


t 

J e-12(t-8) y (s) dv(s) +0 ast—+ oo 

t 

0 
follows by using the same arguments as in Raghavendra and Rao‘ by using property 
(Gi). We next show that 


t 
f e-ale-*) |h(s) | ds +0 ast oo 
0 


we have 


; t/2 
j e-412(t-8) | h(s) | ds = J eTe80=2) | h (s) | ds 
0 ‘0 


t 

+ f e~4/2(t-8) lA (s) | ds 
t 
0 


t/2 
< | A (0) | e-@/2t J e~9128 ds + h (t/2) | ; e~@!2(t-8) | ds 
t/2 


from which it is easy to see that the ri : 
ght hand side of th : 
ZEIO as f > ©0, ¢ above inequality goes to 
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Hence 
| y(t, to, Xo) | > O ast > ©. 


This also proves that if g (ft, 0) = 0 and h(t) = 0, then the null solution of (4.1) is 
eventually uniformly stable, which completes the proof. 


5. A New RESULT 


In the present section we give a new result on necessary and sufficient conditions 
on stability of a non-linear differential system with impulsive perturbation. 


Theorem 5.1—Let x = 0 be uniformly asymptotically stable. Let the condition 
(i) be satisfied for f (t, x) as in Theorem 3.1. Also let g(t, x) satisfy the condition 
(Gi). Then for the following system 


Gx = f(t, x) + g(t, x) Dv + h(t) pt Set) 
the conclusions of Theorem 3.1 hold if and only if h (t) satisfies condition (Go). 


Proor : Let A satisfy condition (G2), then the proof of the Theorem 3.1 
establishes this one way conclusion. 


Conversely, suppose / does not satisfy the condition (G2). Then there exist 
7 > 0 and sequences {tn} and {@n} with 0 < gn < 1 for every n and tn — 0° as 
h — co, such that 


t,*@ 


n 
ie dt | > 7 sd a 
tn 
for every n. Suppose that there exist some f9 > 0 and some Xo for which the solution 
y (t, to, Xo) of (5.1) satisfies | y (ft, to, Xo) | > 0 ast > 9. Choose T so large that 
t > T implies 


| y(t, to, x0) | <4 ({z}-? 


where L is the Lipschitz constant as in condition (!) of Theorem 3.1, we may assume 
without loss of generality that L > 1. Also from condition (G1) it follows, choosing 
n large enough that we can get 


tytn 


[ rOaM< | 


tn 


choose n so large that tn > 7, then we have 


t, +0 tnt®n 


1p” n(ayatl =f F(t 7 Ct tos x0)) dt 


‘n 


(equation continued on p. 454) 
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t+e 


4" g(t,y Ct toy x0) dv (2) 
tn 
tyten 
_ i) y’ (t, to, Xo) dt | 
tytn tytn 
< J L | y(t, to, xo) | dt + J y (t) dv (t) 
n n 


+ | y(tn + On, tn, xo) | + | y (tn, to, Xo) | 
< 9/4 + 7/4 4+ 9/4 + 9/4 = 0 


which contradicts condition (5.2) hence | » (t, to, xo) | can not tend to zero as tf > oo, 
Completing the proof of the theorem. 
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The propose of the present paper is to find the relation between y-curvature 
tensors with respect to CI ofthe Fins!er /spaces (M™, L) and (M”, L*) 
where L* (x, y) is obtained \rom L (x, y) by the transformation. 
\ ; A t ° 
L*? (x, y) = L* (x, p) HU Y*)* 


where X; (x, y) is an A-vector in (M”) L). The relation in n-fundamental 
forms of tangent Riemannian hypersurfaces of (M™", L) and (M", L*) have 
also been obtained. } 


1. INTRODUCTION 


Let Fm = (M”, L) be an in n-dimensional Finsler space, where M” is an n-dimen- 
sional differentiable manifold and L (x, y) is the Finsler fundamental function. 
Matsumoto’? introduced the transformations of Finsler metric : 


L (x,y) = L(x, y) + My tii) 

L*2 (x, y) = L2 (x, y) + (% y')? at lee) 
and obtained the relation between the imbedding class numbers of tangent Riemannian 
spaces to (M”, L), (M™ L) and (M”, L*). It has been assumed that the functions X% 
in (1.1) and (1.2) are functions of coordinate only. Since a concurrent vector field 
is a function of coordinate only, assuming Xi as a concurrent vector field, 
Matsumoto’ has studied the R3-likeness of Finsler spaces (M*, L) and (M”, PAR 
Singh and Prasad!4 and Prasad et al.1° generalized the concept of concurrent vector 
field and introduced the semi parallel and concircular vector fieldswhich are functions 
of coordinate only. Assuming X% as a concircular vector field, Prasad, Singh and 
Singh!9 has studicd the Ri-likeness of (M”, L) and(M”, L). 


If L (x, y) is a metric function of Riemannian space, than L (x, y) reduces to the 
metric function of Rander’s space. Such a Finsler metric was first introduced by 
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Rander’s!2 from the stand point of general theory of relativity and applied to 
the theory of electron microscope by ingraden? who first named it a Rander’s 
space. The geometrical properties of the space have been studied by various authors. 
Numata® studied the properties of (M”, L) which is obtained from Minkowskian space 
(M”, L) by the transformation (1.1). In all these works the functions Xi are assumed 
to be a function of coordinate only. 


Izumi! while studying the conf¢! 4 a Jormation of Finsler spaces, introduced 

the A-vector Xi which is v-covarian{ ... ”.h respect to Cartan’s connection C T° 
c h hab; : ; 

and satisfies L Ci, Xn = p hij. Thus th- “-vector Xi is not only a function of coordi- 


nate but it is a function of direction/4 = nent satisfying L 01 Xi = P hij, Prasad 
etal.\1 has obtained the relaticn | y (¢, 469. imbedding class numbers of 
(M", L) and(M*, L) where L (x, y) d from L (x, y) by the transformation 
(1.1) under the assumption that X% is an in (M”, L). 
2. THE FINsLeR Space ,(M”, LS} 
Let Xi be a vector field in the Finsler space (M, L). If Xi satisfy the condi- 
tions of 


Xely <0 t (2.1) 
je Ci, Xn = P hij »-. (2.2) 


then the vector neld Xi is called an /A-vector!. Here | ; denote the v-covariant differen- 
bs 4 : : a 
tiation with respect to Cartan connection C I’, C, ; is the Cartan’s C-tensor, hij is the 


angular metric tensor and P is function given by 


] 
= ee i 
p a i Guxs ee pty 
where C' is torsion vector Cjz g/*. The first of the following lemmas has been proved 
in! while the other two is a direction consequence of the definition of h-vector. 


E 2.1— i : i 
emma 2.1—If Xi is an A-vector than the function P and X? = Xi — pls are 


independent functions of y. 


Lemma 2.2—The magnitude X of an h-vector Xi is independent function of y- 


Lemma 2.3—For an h-vector Xi we have SnijkX = O where Shijk is v-curvature 
tensor of Cartan’s connection CT. 


Let Xi is an A-vector in the Fnsler space (M", L) and (M”, L*) be another Finsler 
space whose fundamental metric function L* (x, y) is defined by 


L*2 (x, y) = L? (x, y) + B2 (x, y) ...(2.4) 
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where B (x, y) = Xz y!. Since Xi is h-vector from (2.1) and (2.2) we get 
8; Xi = L-1p hus 


which after using the indicatory property of hi; a; 6 = Xj. Thus differentiation of 
(2.4) with respect to )* gives 


L*[* = Lk + BX sy 


where /* = @; L* is the normalized element of support in(M”, L*). The quantities 
of (M”, L*) will be denoted by star letters 
Since 0; 1j = L-1 hij, differentiation of (2.5) with y/ and application of (2.4) 


give 


hr. ras bee is =ohiy + hk lj+Xi Xj ..-(2.6) 


where 


iu (1 ae Er). | atari) 
Hence we have 


vig =gii t (1 -—a)h ly + Xt Xj. He 48.9 | 


From (2.8), the relation between the contravariant components of the fundamental 
metric tensors can be derived as follows : 


(1 — so) (X2 + a) 


gr = ogi OOP i ys pv xy + : Ip 
+ <x! xd ...(2.9) 
where 
A= {So * a3 ¢ Ba } (2.10) 
and X is the magnitude of the vector X‘(=g‘/ Xj). 
From the Lemma 2.1 and relation (2.7) we get 
dio = mi Ree Ly 
where 
pee ye En, sah, 12) 


B. N. PRASAD AND LALJI SRIVASTAVA 


458 
Since 0; kk = L-1 Aus, differentiating (2.8) with respect to y* and using (2.4), (2.11) 
and (2.12) we get 
ee + hyk mi). es sy 
Ch eT Cte + 5 (hts me hrs mj + hye mi) 
From the definition of mi, it is evident that 
(a) mlé =0 
B2 
(b) mt Xt = mm! = X2 — 72 
(c) hig mt = hij X* = mj 
...(2.14) 


(a) Ci, mn = LP hy. 
From (2.9), (2.2), (2.13) and (2.14) we get 


ehh. Sh p a “ 
CH= Cyt spent +H me + hy 





+ Bl fete ( 2 z)jan enn 
taste 





We shall now find the v-corvature tensor Sa je «(Of (M*, L*) which with respect to the 


Cartan’s connection C I’, is defined as 


m ™ 
Sie Chm Cty — Cay, Cte - (2.16) 


Firstly from (2.13) and (2.15) we have 


Cr Ce = Chkm Ci + ay hig hnk + i (Cijk mn 


+ Ctjh mk + Cink my + Cink mi) 


p26 p2 p2 
+ hnk mi mj | aes 2126 + 2L2 Xr 


R2 yt p2 p2 
x2— = ee er ee 
x ( 72 + hij mn mk site + TET ° 


p2 
+4(x2- rz )t|+ Tit (hin mi mk + hin mj 


(equation continued on p. 459) 
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x mk + hig mi mn + hik mj mn) + mi mj Mh Mk 


...(2.17) 


p2 
L2A 


where 








et Ghia al el Fal A 1S ea Pa fo Soe Be 
Paras A RESTS ene area eh cs(* - 4) 


ee BP \. 
ae (xe — E je. (248) 


Thus from (2.16) we get 


Sie = 6 Snizke + hij dnk + hak + hank diz — bik dng — hnj diz 
..(2.19) 
where 
dij = } a1 hts + a mi mj ...(2.20) 
ao ey PS (e-E)-a On 
3. HypersSuRFACE OF (M2, L) 

Let (M"-1, L) be a hypersurface of (M”, L) given by the equation 
x? = x! (u*). oaaloale 
Let us suppose that the functions (3.1) are atleast of class C? in u* and the pro- 
jection factors a = es are such that their matrix has maximal rank n — 1. The 


fundamental metric function L (, v) of the hypersurface is given by 
L (u®, v®) = L (xt (u*) BY v*) 
where v* is the element of support for the hypersurface for which 
yo = BL. 
Thus if /* denote the normalized vector using the element of support then 


Riad Be [*, soe(3.2) 


If gnj (x, y) denotes the metric tensor of (M”, L), the induced metric tensor of (M"-1, L) 
is given by 


h >of 
gap (u,v) = gnj (x,y) B, Ba. ..(3.3) 
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The inverse of (3.3) is denoted by g%8 (u, v) by means of which we define the quantities 
j 
BY (n, v) = g*®(u, v) ges (x, y) Bg - (3.4) 


The unit normal vector N/ (x, y) of (4-1, L) is determined by the relations 
ghj (x, y) B* BNI (x, y) = 0, gnj (x, y) N* (x, y) Ni (x, y) = 1. ...3.5) 
We have the following identity from (3.3), (3.4) and (3.5) 


x pl j & pe sanias0 
By Bz, = 83, B, By, + NIN =8, (3.6) 


where 

Ni = gij (x, y) NI. 
If Cnjz (u, y) denotes the (A) Ay-torsion tensor of (M”, L) the induced (h) hAy-torsion 
tensor Cupy (u, v) of (M”-1, L) is given by 

Copy (u, ¥) = Crs (x,y) B, By Be aaa 
from which we obtain 


a -s « j j 8 
Coe BC BR, 3.8) 


The relative v-covarient derivative of the projection factor B; with respect to the 
induced Cartan connection / C I is defined as® 


é t t t h k 
Bye titara Bing Cy yt Coat OB Ae, ...(3.9) 


This tensor is normal to(M"-1, L). Therefore we may write 


t . 
Bay — May Nt. ...(3.10) 


From (3.9), it is clear that Mpy is symmetric in 8 and y and it may be written as 


Mpy = Cex Nt By BY Re SG) 


The tangent vector-space Moon to M”~1 at every point x! (= u*) of the hypersurface 


is considered as the Riemannian space (M* ‘ » §z) with the Riemannian metric Qc = 


8x8 (u,v) dv* dvB. The components of the (h) hy-torsion tensor Cay will be the Chri- 


stoffel symbols associated with gz. If M’, is the tengent vector space to Mat xt 
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(= u*) the (M* : gr) will be the hypersurface of (M* , Zz) given by (3.2a) where 


22 = gii (x, y) dy! dy9 is the Riemannian metric of M* The quantities Mgy given in 
(3.11) will be considered as the coefficients of second fundamental forms of tangent 


Riemannian space (M ? a ga). 


In general the coefficients of the rth fundamental forms of (M*-1, gz) are defined 


asi8 
C(1)a8 = SaB 
C(2)08 = Map 
Cap = Claes Mg (2<r <2) 
where 
Mg — 8** Mas. 
4. h-Vector Fieips 1n (M"-1 L) 
At the point of (M”-1, L), the vector field Xi may be written as 
Xe Xen B, ip Ne (4.1) 
where 
(a) Yoo XB (by p= XN (4.2) 


Since Xile = Xtlj By, we have from (4.2 (a)) and (3.10) 


Yale = Xtly BL By + » Mas. .»(4.3) 

From (3.1!) and (3.6) we get 
. . Ve aeons aa .. (4.4 
DC yy tem LC) Xs By BL — Lv Mav. (4.4) 


If Xi is a concurrent vector field in Fn then in view of (2.1), (2.2), (3.2b) and (3.5) 
eqns. (4.3) and (4.4) reduce to 


XalpB = b M.? L Coy Xe. = Phay — Lp May 


These relations yield the 
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Theorem 4.1—If Xi is an h-vector field in (M", L) the vector field X. = Xt 
B‘ is also an h-vector field in (M"-1 L) if and only if 


(i) Xi is tangential to the hypersurface (M"~1, L) 


or 
(ii) Mas = 0. 

The hyperplane of first, second and third kinds are defined®. Ina hyper plane of 
third kind Mz vanishes®. Thus : 


Theorem 4.2—If Xi is an h-vector field in (M”, L) then vector field X% Bt is also 
an hA-vector field in a hyperplane of third kind. 
In the following we assume that Xi is tangential to (M@”~!: L), so that 


(a) Xi = Xa BX (b) Xi = X* Bt (4.5) 


a 


where 
X* = g*B Xp. 


5. THE n-FUNDAMENTAL Forms OF HYPERSURPACE OF (M®, L*) 


Let (M""!, L*) be a hypersurface of (M™, L*) given by the same equation (3.1). 
The relations (2.8), (2.9), (3.2b), (3.3) and (4.5) yield 


Seg = 9 8ab(1 — 6) le Ip + Xo Xp (5.1) 
tap — 1 g®B __ igh: a VB B o 
g ae LA (i= X8 4 [8B ye) 
(X2 + o2) l 
i ewe Ls sf 8 | Pa 
+ ( s) 7 I= JB 4. x= x8. aula e) 
From (3.2b), (3.3) and (2.5) we also have 
; 
has = hey BL Ba , hes Bi Ni=0 is, (5,3) 
and 
Map = Xy — B/L l, = me Bt ; ...(5.4) 


It is to be noted that if N* j i 
1S a unit ies it i 
Aleta normal vector to (M*~1, L*) then it is not normal 


We may write 
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Nti = N2 B. 4+ DNi. (5.5) 


~ 


To obtain N® and D we use (3.8), (3.5) and (5.1). Thus we get 
g*.e N* + Du Xa = 0 (5.6) 
get N= NB + D2 0 + » (2D Xa N* + D* un) = 1. 6.7) 
If Xi is tangential to (M”-1, L) then »p = 0 and eqns. (5.6) and (5.7) give N* 
= 0, D = og 1/2 
From (5.5) it follows that 
N*t = o-li2 Ni, ...(5.8) 
Hence we have the following : 


Theorem 5.1—Let (M, L*) bea Finsler space obtained from (M™, L) by the 
transformation (2.4). If (M?~1, L*) and (M-1, L) are hypersurfaces of these spaces 
and Xi is tangential to the hypersurface (M*-1,L) then the vector normal to (M”-1, L) 
is also normal to (M"~1, L*). 


Now we estoblish the following : 


Theorem 5.2—Let (M™, L*) be a locally Minkowskian space obtained from lo- 
cally Minkowskian space (M”, L) by the transformation (2.4). Let (Mn-1, L*) and 
(M-1, L) be hyper surfaces of (M™, L*) and (M’, L) resspectively. If X% is tangential 


to the hypersurfaces (Mn-1, L) and (M* , &x), (M*, a ), (Mm, , 22), (mM, : g° ) 
are tangent Riemannian spaces to (M, L); (M", L*), (Mn-1, L), (Ma? L*) respec- 
tively, then we have the following : 


(i) Second fundamental forms of (me, gz) and (M" ie g*) are proportional 
(ii) Every asymptotic direction of (m,-* _ gz) is asymptotic direction of (mM, * 
94 


(iii) The rth fundamental tensors of (M,* , Zax) and (Mm, i g* ) are related 


by 
r-1 
C* (ap = 08 7!? [Cirae + S Pye O(rem_a)a],3 <r Qn (5.9) 
m=2 
where 


ce] 
- — ee 2 ee re eo | .. (5.10) 
ye ie C(m)as 
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Rom = [Poems X85, 2<Sm<n—1 (5.11) 

O(2)0 = P(2)a ...(5.12a) 
f—1 

O(rn = Pir)e + 2 Rim) Q(rgi-m)x, 3S r<cn— | oy fH PA) | 
mM=Z 


Proor: (i) If Xt is tongential to the hypersurface to (M”-1, L) then » = 0 
and hecce 


m Nt = 0. 
Thus from (2.13, (3.10), (5.3), (5.7) it follows that 
Mpy = 31/2 Mpy. 5.19) 
This proves (i). 


(ii) A direction ** forwhich Map t* t® = 0 is said to be an asymptotic dire- 
ction. In view of this dafinition and (5.13) we get (ii). 


(ii) The validity of relation (5.9) is established by induction. 


Since Cijx is an indicatory tensor from (3.2b) and (3.11) it follows that Mpy 
1 = 0. Hence from (3.12) we get 


Ciray 1¥ = 0 = Ciryay 1? 2 y En. ...(5.14) 
Hence from (5.10), (5.11), (5.12) we get 


P(r) I* = 0, Oye =O, 2Er<n, (5.15) 
From (5.2), (5.13) and (5.4), we get 


oa. ~ ia ss, ma (l — o)B a J@ o 
Mayne ag [ Mm a | May XY. ...(5.16) 
The relations (3.12), (5.10), (5.14), (5.15) and (5.16) yield 
o 
Cae = Coan + ,f ) Pore Pea. (6.17) 


From (5.12a) and (5.17), it is evident that (5.9) holds for r = 3. 


Fora given fixed 
value of the integer s with3 <s <n — 1, we have 


C* (841) 0R a ee M;° ; 


++0(5.18) 


Now let us suppose that (5.9) is valid for s = 34 oe 


.. fr, SO thai we ca ; 
in the form N write (5.18) 
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$_} 
CCet1)a8 = 9 Cloyan + SS P(m)p Q(s41-m) ] M;* 


m=2 


which in view of (5.10), (5.11), (5.12a), (5.14), (5.15) and (5.16) gives 


s-1 
Coet)ab = a(2-M Coo a)as aE > P(m41)8 Q(s41_m)« 
m=2 


8-1 
+ P(2)p {P(2)s + ee) Q(s,1-m)«} ] 
m= 


& 
= a(2-8)I2 [C(syo)an + % peed y QO(s,2-m)a]. 
iol 


This show that (5.9) is valid for r = s + |, which completes is proof of (iii). 
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REMARKS ON SUBMANIFOLDS OF CODIMENSION 2 OF AN EVEN- 
DIMENSIONAL EUCLIDEAN SPACE 
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Department of Mathematics, Faculty of Science, Hiroshima University 
Hiroshima 730, Japan 


(Received 10 May 1988) 


We determine submanifolds M of codimension 2 in an even-dimensional 
Euclidean space E in relation to the integrability of the almost complex 
structure J on M X R! X Rl. 


INTRODUCTION 


Blair et al.1 introduced the so-called (f, g, u, v, A)-structure which is naturally 
defined in a submanifold of codimension 2 of an almost complex manifold. An even- 
dimensional sphere of codimension 2 of an even-dimensional Euclidean space is a 
typical example of a manifold admitting the structure. Okumura2, Yano and Okumura’ 
studied submanifolds admitting a normal (J, g, u, v, A)-structure. 


In the present paper, we shall study submanifolds of codimension 2 in an even- 
dimensional Euclidean spece in relaticn to the Nijenhuis tensors formed by the in- 
duced (f, g, u, v, A)-structure on the submanifolds. 


1. SUBMANIFOLDS OF CODIMENSION 2 OF AN EVEN-DIMENSIONAL 
EUCLIDEAN SPACE 


Let E be a (2n + 2)-dimensional Euclidean space and X the position vector of 

a point of E from the origin of E. Since E is even-dimensional, it is regarded as a flat 
Kaeblerian manifold, that is, there exists a (1, 1)-tensor field F satisfying 
F2=—J], FY-FZ=Y.2Z Ft) 


for any vectors Y and Z and 


VF =0 rea bee) 
where J denotes the identity transformation, a dot the inner product in £ and V the 
Riemannian connection of E. 


Let M bea 2n-dimensional submanifold of EF covered by a system of local co- 


ordinate neighbourhoods {U ; x*}, where and in the sequel the indices h, i, j, k 


aap 
over the range {1, 2, ..., 2n}. 
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By putting Xi = ¢@ X (0% = 4/dx'), Xi are 2n linearly independent local vector 
fields tangent to M and the induced Riemannian metric of M is given by gy = X;. Xi. 


We denote by C and D two mutually orthogonal unit normals to M such that Xt, 
C and D form the positive orientation of E. 


The transforms FXi, FC and FD of Xi, C and D by F can be expressed as 


FXi =f, Xn, +uC+wD ae =o) 
FC =-— ut Xx + AD 
Ki DisenvtXe —.C (1.4) 


where fe are components of a tensor field of type (1, 1), w and vw 1-forms, A a func- 


tion of M, u? = ui gi* and v’=% git. 


Applying F to (1.3), (1.4) and (1.5) and taking account of (1.1) and (1.3)~ (1.5) 
we have2’4. 


rs 74 =— 3, + usu® + vi vi ...(1.6) 
ut f =A vi, vif; = — Au =i) 
uui=1—=—wWv', uw vi = 0 cag F,-}) 
i fy gkh = gt — usu — vj ..-(1.9) 


h ; vey : 
and fii = f, ght is skew symmetric in J and j 


The totality (/, g, u, , A) satisfying eqns. (1.6) ~ (1 .9) is called an (f, g, u, v, A)- 
structure!. 


The equations of Gauss and Weingarten are given by 


wi Xi = ha C+kn D ..-(1.10) 
eer (lit) 

vi C=—h, Xn+ UD (1.11) 
s (1.12 

9 D=— kX -— GD (1.12) 


where hy: and kj are the second fundamental tensors of M with respect to the normals 
C and D respectively, and /j the third fundamental tensor of M in E and y the Rieman- 
nian connection of M. 
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Differentiating eqns. (1.3) ~ (1.5) and taking account of (1.2), (1.10), (1.11) and 
(1.12), we have?’4. 


Vi fa = — hau? + hi ii Ky ky vi (1213) 
Viui = — Air ts —Aky + i vi .--(1.14) 
vin = — kik fp + Ahi — Lue ruts) 
Wid = — Ay ve t+ ky ue. ...(1.16) 


2. NUENHUIS TENSORS 


We consider the product manifold M x R1! x R1, R! being a 1-dimensional 
Euclidean space, and we denote it by M. The indices A, B, C, ... will run over the 
ranges 1, 2, ...,27, 2n + 1, 2n + 2, and * stands for 2m + 1 and + for 2n + 2. If 
we define on M a tensor field J of type (1,1) by local components 


[ if ui VE 
A | | 
(Jp) heard ay ee £2 
| 
| —v — AD | 
in each {U x R! x RI, x4}, then we can see that J? = — Jon M, that is, {M, J} 


becomes an almost complex manifold. The components of Nijenhuis tensors of the 
almost complex structure J are given by 


h k A k h k BA ua 
Nn = t, Us, — Sy Ob fi Cie cee ee sol enh 
+ (05 ut — Of uj) uh + (8; ve — % vj) vA 


* _ gt k BLS ae 
Ni =; Okut — f, Ok uj — (0; ie OU Ss ) uk 


+ A (0; vi — O% 5) (223) 
+ k k 
Nii = f; Oevi — f, Ok vi — (45 fi — & fi) ve sae( 24) 
— Xr (j ut — 8 uj) 
h k 
Nig = — f, avd — ve ORS? + (Os uk) f? + (ea) vr ...(2.5) 


h k h 
Nige ~~ Sy Oe + ve ORS; + (es vi) FE — (sa) 


h 
N = yk 1s, 
op 6 Ne OM ME eee (UR 
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Ne = —u¥ ev —f, aA — (Uk) ve -..(2.8) 
Ne, = — uk de us — (& uk) ux — A (4 A) wih 29) 
Nya = — fy (er) — vE (8x uj) + (2s VE) ue ...(2.10) 
Nie = vk (ae vj) + (8; VF) ve + A (2 A) me 2) 8 
Nive = — ut (eA) (2.12) 
Nave =k (28a). (2.13) 


The Nijenhuis tensor of an almost complex structure J satisfies the identity® 
7A E Byd> & 
Noe J ae Nes Jy =. 


Hence, substituting (2.1) into (2.4) and taking account of the algebraic properties of 
(f, g, u,v, A)-structure, we can obtain the 

Proposition 2.1—Let there be given an (f, g, u,v, A) structure on M and let the 
function A (1 — A2) does not vanish almost everywhere. If Ni, A. re and N+ vanish 


on M, so do the other components of the Nijenhuis tensors of J. 


3. SUBMANIFOLDS WITH VANISHING NUJENHUIS TENSORS 


Let M be a submanifold of codimension 2 in a (2n + 2)-dimensional Euclidean 
space E. Then M admits an (f, g, u, ¥, A)-structure and the equations of Gauss. 
Codazzi and Ricci on M are given by 


Kr, = hy hie — i hee + kn Ky — ky kext ..B.1) 
vehi — Vile — Ikkan + ly kit = 0 ...(3.2) 
ve kis — 07 ket + le hn — yp hit = 0 natso) 
vile — els + Hy kw — hi ku — hy ky = 0 (3.4) 


respectively. 


Assume that the function A(1 — A2) does not vanish almost everywhere on M. 
In the case where theconnection induced in the bundle of M has null curvature, we 
say that the normal connection is trivial. 


The following lemma is well known4. 


470 BYUNG HAK KIM 


h . i : 
Lemma 3.1—If the tensor N it vanishes on M and the connection induced in the 


bundle is trivial, then there exist two constants « and B such that 


, ae ; 
h, ul = at, h, vi = ayt 


hut 
J 
k, u 


l 


Bui, kK, vi = B yt 
h, hih = « hji, ky kin = Bk 


and hi and ke commute with ie 
On the other hand, the tensor fields (2.2) ~ (2.4) can be rewritten as 
h t ,h t ,h h ,h h 
Ny CF; Ay — hy fy De — (fp yok, hy Sy Dus 


h 


+(fik, —KSf, u—-Uike —k fh yy 
+ (lj vi — he vj) uh — (jy ut — ik uj) vv’. 


Ni = le (f; Vi — fe yy +A(k ui — I; ut) 


k k k 
— (hast k, aii, uj — ke vj) ue 


Ne =e (ft uy — f* wt) — Ab — bevy) 


— VE (h, ui + ki “ — he uj — ke vj). 


If os = 0 in addition to the assumptions of Lemma 3.1, then we have 
Buy vi — Bu vy = 0. 
Similarly, Ne = 0 implies 


ovum —avwuy = 0. 


(3.5) 
...(3.6) 


ato) 


2435) 


(3.9) 


med fe 2 (8) 


PoC G4 


...(3.12) 


h 
Lemma 3,2—If the tensor N,, vanishes on M and the connection induced in the 


normal bundle is trivial, then the 


(1) The components Ni vanishes. 


(2) 8 =0. 
(3) kw = 0. 


following condition are equivalent to one another : 
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In this case, the equations of Gauss and Weingarten reduce to 
wi Xi = hat Ci wi C= - hy Xn. Yi D= 0 


respectively, and consequently D is a constant vector and we have 
vi(X - D)=0 
that is, X - D = constant. Thus we have 


Theorem 3.3—Let M be a submanifold of condimension 2 of E2"*? and suppose 
the connection induced in the normal bundle of M is trivial. If Ni, and NS vanish, then 


M isa hypersurface of a hyperplane E2"*1. 
We prepare the following : 


Theorem A4—Let M be a 2n-dimensional complete differentiable hypersurface in 
a (2n + 1)—dimensional euclidean space E’. If the component Ni, vanishes on M, 


then M is a product of a sphere and a plane. 
Combining Theorems 3.3 and A, we obtain: 


Theorem 3.4—Let M be a complete submanifold of codimensian 2 of E2n*2 such 
that the connection induced in the normal bundle of M is trivial. If Ni, and N i va- 


nish, then M is a product of a sphere and a plane. 


Assuming Ne 0 instead of Ape = 0, we can obtain the same conclusion as 


h 
Theorem 3.4. Therefore if the components J, , Nas and Ne 
the connection induced in the normal bundle is trivial, then we get 


vanish identically and 


yw Xi = 0, C= 0, D=0 
and 
Theorem 3.5—Let M be a submanifold of codimension 2 of E2+2 with trivial 
normal connection. If the Nijenhuis tensors Ni, , ‘Ny and N fi vanish on M, tben 


M is a plane of codimension 2. 
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A function f: X > Y is said to be almost irresolute'? if foreach x € X and 
each semi-neighbourhood V of f(x), the semi-closure of f-1(V) is a semi- 
neighbourhood of x. In this paper, we obtain several characterizations of 
almost irresolute functions and investigate the relationship between such 
functions and some weak forms of irresolute functions. We also improve on 
some results established by Dube et al.°. 


1. INTRODUCTION 


Crossley and Hildebrand® defined irresolute functions by utilizing semi-open sets 
due to Levine!?, Recently, as weak forms of irresoluteness, weak irresoluteness9, 
6-irresoluteness9, almost irresoluteness?! and quasi irresoluteness® have been defined 
and investigated independently. However, it will turn out that these four weak forms 
of irresoluteness are equivalent. On the other hand, Dube et a/.19 have introduced 
the notion of almost irresolute functions which is independent of that of almost 
irresolute functions in the sense of Thakur and Paik?!, 


The purpose of the present paper is to investigate almost irresolute functions in 
the sense of Dube et al.1°. Note that, after section 3, ‘almost irresolute” always means 
“almost irresolute” in the sense of Dube etal. In section 2, we point out that 
weak irresoluteness, @-irresoluteness, almost irresoluteness in the sense of Thakur 
and Paik and quasi irresoluteness are all equivalent. In section 3, we present several 
characterizations of almost irresolute functions in the sense of Dube etal. In 
section 4, we investigate the relationship among semi-continuity, quasi irresoluteness 
and almost irresoluteness in the sense of Dube et al.1°. In section 5, we introduce 


* This research was supported by a grant from the C.N,A. (G.N.S.A.G.A.) and M.P.I, through 
**Fondi 40%”’. 
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and investigate semi preopera functions. The last section concerns with strongly 


semi-closed graphs and contains several improvements on results established by Dube 
et al.9. 


2. PRELIMINARIES 


In this section, we will point out that almost irresoluteness?!, weak irresolute- 
ness?, 9-irresoluteness? and quasi irresoluteness® are all equivalent. 


Throughout the present paper, (X, +) and (¥,<¢) (or simply X and Y) denote 
topological spaces on which no separation axioms are assumed unless explicitly stated. 
A subset S of X is said to be semi-open!? if there exists an open set U of X such that 
UCSC CI(U), or equivalently if S C Cl (Int (S)), where Cl (S) and Int (S) denote 
the closure of S and the interior of S, respectively. A subset S is called a semi-neigh- 
bourhood? of a point x of X if there exists a semi-open set U such thatx € UC S. 
The complement of a semi-open set is called semi-closed. A semi-closed and semi- 
open set is said to be semi-clopen. Fora subset S of X, the intersection of all semi- 
closed sets containing S is called the semi-closure of S (Crossley and Hildebrand®) and 
is denoted by s Cl (S). The semi-interior of S, denoted by s Int (S), is defined by the 
union of all semi-open sets contained in S. The family of all semi-open sets of X is 
denoted by SO(X). For each x € X, the family of all semi-open sets containing x is 
denoted by SO (X, x). A subset S is said to be regular-open if S = Int (C1(S)). A 
subset S is said to be regular semi-open? if there exists a regular open set U of X such 
thatUC SC CI(U). 


The following two lemmas are due to Di Maio and Noiri’. 

Lemma 2.1—The following are equivalent for a subset A of a space > 
(a) A is regular semi-open. 

(b) A = s Int (s Cl (A)). 

(c) A is semi-clopen. 

Tenma?2—If-A © SO (X), thens Cl (A) ts semi-clopen. 


Definition 2.3— A function Pesta Ate oe tae «) is said to be irresolute® (resp. 
semi-continuous!2) if f-1 (V) € SO (X, «) for every V € SO(Y, cs) (resp. V € 9). 


Definition 2.4— A function f: X > Y is said to be (a) almost irresolute”* if 
-1(V) € SO (X) for every regular semi-open set V of Y; (b) weakly irresolute® (resp. 
§-irresolute®) if for each x € X and each semi-neighbourhood V of f (x), there exists a 
semi-neighbourhood U of x such that f(U) C s Cl (V) (resp. f (s Cl (U)) eae 9 81 4 8) » 
(c) quasi irresolute® if for x € X and each V € SO (Y, f(x)), there exists UE SO(X, x) 
such that f(U) C s Cl (V). 


Lemma 2.58—The following are equivalent for a function /f: ¥ > Y: 
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(a) f is quasi irresolute. 


(b) For each x € X¥ and each V € SO (Y, f (x)), there exists UE SO (X, x) 
such that f (s Cl (U)) C s Cl (V). 


(c) f-1 (V) is semi-clopen in Y¥ for every semi-clopen set V of Y. 

(d) f-1(V) C s Int (f 71 (s Cl (V))) for every V € SO (¥). 

(e) sCl(f-1 (V)) C f-1(s CI(V)) for every V € SO (Y). 

Theorem 2.6— The following are equivalent for a function f: X¥ -> Y: 
(a) f is quasi irresolute. 

(b) fis weakly irresolute. 

(c) fis 6-irresolute. 

(d) f is almost irresolute. 

Proor: This follows from Definition 2.4, Lemmas 2.1 and 2.5. 


Remark 2.7 : It is shown in Di Maio and Noiri® that semi-continuity and quasi 
irresoluteness are independent of each other and they are implied by irresoluteness. 
3. CHARACTERIZATIONS 


In this section, we obtain several characterizations of almost irresolute functions 
in the sense of Dube ef al.1° A subset S of a space X is said to be preopen!5 if 
S Cc Int (Cl (S)). 


Definition 3.1— A subset S of aspace X is said to be semi-preopen! if there 
exists a preopen set U in X¥ such that UC S C Cl (UV). 


The family of all semi-preopen sets in X is denoted by SPO (X). The comple- 
ment of a semi-preopen set is called semi-preclosed}. 


Lemma 3.21—The following are equivalent for a subset A of a space X, 
(a) A € SPO (X). 

(b) A C Cl (Int (C1 (A))). 

(c) A C s Int (s Cl (A)). 

Proor: This follows from Theorems 2.4 and 3.21 of Andrijevic!. 


Definition 3.3— A function f: ¥ > Y is said to be almost irresolute!® if for 


each x € X and each semi-neighbourhood V of f (x), s Cl (f-1(V)) is a semi-neigh- 
bourhood of x. 


Henceforth, ‘‘almost irresolute” always means “almost irresolute” in the sense 
of Definition 3.3, that is, Dube et a].10, 
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Theorem 3.4— The following are equivalent for a function f: X > Y: 
(a) fis almost irresolute. 

(b) f-1(V) C s Int (s Cl (f-1 (V))) for every V € SO (Y). 

(c) f-1 (VV) C Cl (Int (Cl (f 1 (V)))) for every V € SO (¥). 

(d) f-1(V) € SPO (X) for every V € SO (Y). 


Proor : (a) = (b): Let V © SO(Y) and x € f-! (V). Since V isa semi-neigh- 
bourhood of f(x), s Cl (f-! (V)) is a semi-neighbourhood of x and hence there exists 
U € SO (X, x) such that UC sCl(f-1(V)). Therefore, we have x € UC sInt 
(s Cl (f-1(V))). This implies that f-1 (V) C s Int (s Cl (f-1(V))). 


(b) > (a): Let x € X and V be any semi-neighbourhood of f(x). There exists 
W € SO(Y, f (x)) contained in V. Therefore, we obtain 


x € f-1(W) C s Int (s C1 (f-1 (W))) C sC1(f-4 (W)) Cs Cl(f4 (V)). 
This implies that s Cl ( f-1 (V)) is a semi-neighbourhood of x. 
It follows from Lemma 3.2 that (b), (c) and (d) are all equivalent. 


Theorem 3.5— A function f: X > Y is almost irresolute if and only if 
f(s C1 (U)) C s C1(f (U)) for every U € SO (X). 


Proor: Necessity— Let U € SO(X). Suppose that y & s Cl ( f(U)). There 
exists V € SO(Y, y) such that VM f (U) = ¢; hence f-1(V) 1 U= 4. Since 
U € SO(X), we have s Int (sCl (f-1 (V)) AsCl(U) = ¢. By Theorem 3.4, 
fiAWV)AsCi(U)= ¢ and hence V ()\ f (s Cl (U)) = ¢. Therefore, we obtain 
y & f(s Cl (U)). This shows that f(s Cl(U)) C s Cl (f (U)). 


Sufficiency— LetV € SO (Y). Since X¥ — sCli(f- (V)) € SO (X), we have 
f(sci(x—sCl(f1(V))) Csci(f(¥ - scl ( f-1 (V)))) and hence 


X — sInt(sCl(f2(V))) CfA (s Cl (F(X - sCl( fF? (Y)Y)) 
cfi(sci(f(xX —f32M)) chica’ —V))=* — f1(V). 


Therefore, we obtain f-1 (V) C s Int (s Cl(f-? (V))). It follows from Theorem 3.4 
that f is almost irresolute. 


Theorem 3.6—The following are equivalent for a function f: X — ¥. 
(a) f is almost irresolute. 


(b) For each x € X and each V € SO (Y, f (x)), there extsts U € SPO (X ) 
containing x such that f(U) C V. 


(c) f-1 (F) is semi-preclosed in X for every semi-closed set F of Y. 


(d) Int (Cl (Int ( f-1 (B)))) C £74 (s Cl (B)) for every subset B of Y. 
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(e) f (Int (Cl (Int (A)))) C s Cl ( f(A)) for every subset A of X. 


Proor : (a) > (b): Let x € X andV € SO(Y,f(x)). Set U =f-1 (V), then 
by Theorem 3.4 U is a semi-preopen set containing x and f(U) C V. 


(b) > (a): Let VV © SO(Y) and x € f-1(V). There exists U € SPO (X) con- 
taining x such that f(U) C V. By Lemma 3.2, we obtain 


x € UC sInt(sCl(U)) C s Int (s Cl( f-1 (V))) 


and hence f-!(V) C s Int (s Cl(f-! (V))). It follows from Theorem 3.4 that / is 
almost irresolute. 


(a) > (c): This is obvious by Theorem 3.4. 


(c) > (d): Let B be any subset of Y. Since s Cl (B) is semi-closed, f-1 (s Cl (B)) 
is semi-preclosed. By utilizing Lemma 3.2, we have 


X — f-4(s Cl(B)) C Cl (Int (Cl (X — f-1 (s Cl (B))))) 
= X — Int (Cl (Int ( f71 (s Cl (B))))). 
Therefore, we obtain Int (Cl (Int ( f~1 (B)))) C f-1(s C1 (B)). 
(d) = (ec): Let A be any subset of X. We have 


Int (Cl (Int (4))) C Int (Cl (Int ( f~! ( f (A))))) C f-1 (s Cl (Ff (A))). 
Therefore, we obtain f (Int (Cl (Int (4)))) C s Cl (f (A)). 
(e) > (a): LetU € SO(X). Sinces C1(U) = U U Int (C1 (U)) = UU Int 


(Cl (Int (U))), we obtain f(s C1(U)) = f (U) U f (Int (Cl (Int (U)))) C F(U) Us Cl 
(f(U)) = s C1( f(U)). If follows from Theorem 3.5 that fis almost irresolute. 


4. COMPARISONS 
In this section, we investigate the relationship among irresoluteness, quasi 
irresoluteness, almost irresoluteness, semi-continuity and weak quasi continuity. 


Definition 4.1— A function f: X > Y is said to be weakly quasi continuous!9 
if for each x € X, each open set U containing x and each open set V containing f(x), 
there exists an open set G of X such that 6 # G C Uand f(G) C Cl (V). 


Remark 4.2: The following implications hold for properties on a function : 
quasi irresoluteness 
irresoluteness => almost irresoluteness weak quasi continuity- 


NN semi-continuity Sore, 


[t is obvious by Lemma 2.2 and Theorem 3.4 that irresoluteness implies almost 


irresoluteness. The other implications have been shown in section 7 of Di Maio and 
Noiri8. 


ALMOST IRRESOLUTE FUNCTIONS 477 


We shall show that quasi irresoluteness, almost irresoluteness and semi-con- 
tinuity are respectively independent. The following example shows that almost irreso- 
luteness does not imply weak quasi continuity and hence it implies neither quasi 
irresoluteness nor semi-continuity. 


Example 4.3— Let X be the set of real numbers, + the indiscrete topology for 
X and c the discrete topology for X. Let f: (X, t) — (X, o) be the identity function. 
Then fis almost irresolute but it is not weakly quasi continuous. 


Example 4.4— Let X = {a, b, c} and + = {¢, X, {a}, {b}, {a, b}}. Let f: (X, +) 
— (X,+) be a function defined as follows: f (a) = f(b) = aandf(c) =c. Thenf 
is continuous and hence semi-continuous. However, f is not almost irresolute 
because there exists {b, c} € SO (X, t) such that f°! ({b, c}) & SPO (xX, +). 


By Examples 4.3 and 4.4, we observe that almost irresoluteness is independent 
of semi-continuity and also it is independent of weak quasi continuity. The following 
example and Example 4.3 show that almost irresoluteness and quasi irresoluteness are 
independent of each other. 


Example 4.5—Let X = {a, b, c}, + = {¢, X, {a}, {5}, {a, b}} ando = {¢, X, {c}}. 
Let f : (X, t) > (X, a) be the identity function. Then fis quasi irresolute. However, it 
is not almost irresolute because there exists {c}E SO (X , 6) such that f~! ({c}) 
& PSO (X, *). 


It has been shown in Examples 7.2 and 7.3 of Di Maio and Noiri® that quasi 
irresoluteness neither implies semi-continuity nor is implied by semi-continuity. 


5. SEMI-PREOPEN FUNCTIONS 
In this section, we introduce the notion of semi-preopen functions which is 
independent of both notions of preopen functions and semi-open functions. It will be 
shown that every quasi irresolute function is almost irresolute if it has one of the 
following properties : “‘semi-preopen’’, “semi-open” and “‘preopen”’. 


Definition 5.1—A function f: X > Y is said to be semi-preopen if f (U) 
€ SPO (Y) for every U € SO (X). 


Definition 5.2— A function f: X —> Y is said to be semi-open2 (resp. preopen!®) 
if f (U) is semi-open (resp. preopen) in Y for every open set U of X. 


Rose29 called preopen functions almost open and showed thatf: X — Y is 
almost open if and only if f~1(CI(V)) C Cl(f-1(V)) for every open set V of Y. 
Therefore, ‘‘preopen’’ is equivalent to “almost open” in the sense of Wilansky?*. 


We shall show that “‘semi-preopen”, “‘preopen’”’ and “‘semi-open” are respectively 
independent. It follows from Examples 1.8 and 1.9 of Noiri!? that “preopen” and 
“semi-open” are independent of each other. The following example shows that “open 
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does not imply “‘semi-preopen” in general. Therefore, neither preopenness nor semi- 
openness implies semi-preopenness. 


Example 5.3— Let X = {a, b, c, d}, + = {¢, x, {a}, {a, c}} and« = {¢, x, {a}, 
{c}, {b, c}, {a, c}, {a, b, c}, X}. Let f: (X, t) > (X, 4) be the identity function. Then 
fis an open function. However, f is not semi-preopen because there exists {a, 5} 
€ SO (X, +t) such that f ({a, b}) & SPO (X, ¢). 


Let f: (X, t) > (X, ¢) be the function of Example 4.3. Then /~-! is semi-pre- 
open but not semi-open. Moreover, the following example shows that a semi-preopen 
function is not necessarily preopen. 


Example 5.4— Let X = {a, b, c}, + = {4, X, {a}, {a, b}} and « = {¢, X, {a}, {5}, 
{a, b}}. Let f : (X, +) > (X, ¢) be a function defined as follows : f(a) = a, f(b) =, 
and f(c) = 6. Then fis semi-preopen. However, fis not preopen because there exists 
{a, b} € + such that f ({a, 5}) is not preopen. 


Theorem 5.5—The following are equivalent for a function f: ¥ > Y: 
(a) f is semi-preopen. 

(b) f~1 (s C1 (V)) C s Cl( f-1 (V)) for every V € SO (Y). 

(c) f-} (Int (Cl (Int (B)))) C s Cl (f-1 (B)) for every subset B of X. 
(d) f (Int (A)) C Cl (Int (Cl ( f(A)))) for every subset A of X. 

(e) f(U) C Cl (Int (Cl (f (U)))) for every U € SO (X). 

(f) f(U) C s Int (s Cl (f (U))) for every U € SO (X). 


Proor : (a) > (b): LetV € SO (Y) and x & sCl(f-1(V)). There exists 
UE SO(X, x) such that U WM f-1(V)=¢; hence f(U)  V = ¢. Therefore, we 
have s Int(s Cl(f(U))) 1 V = ¢ and hence s Int (sCl(f(U))) Os C1\(V) = ¢. 
Since f is semi-preopen, by Lemma 3.2 we obtain f (UU) NsCl (V) = ¢ and 
U (\ f-1(s C1 (V)) = ¢. Therefore, we have x & f-1 (s Cl (V)) and hence f-1 (s Cl (V)) 
C sCl(f-! (V)). 


(b) > (c): Let B be any subset of Y. We obtain 
f~* (Int (Cl (Int (B)))) = f-1 (s Cl (Int (B))) C s C1 ( f-} (Int (B))) 
C sCi(f-! (B)). 
(c) > (d): Let A be any subset of ¥. Then, we have 
xX — f+ (Cl (Int (Cl ( f (A))))) = f-1 (nt (Cl (Int (Y — f (A))))) 
C sCl(f (Y — f(A))) C sCl(X¥ — A) = X — ss Int (A). 
Therefore, we obtain f(s Int (A)) C Cl (Int (Cl ( £(A)))). 
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(d) > (e): Let U € SO (X). We have f (U) = f (s Int (U)) C Cl (Int 
(Cl ( f (U)))). 


(e) > (f) and (f) > (a): These follow from Lemma 3.2. 


Theorem 5.6— A quasi irresolute function f: X — Y is almost irresolute if it 
satisfies one property of the following : 


(a) semi-preopen, (b) semi-open, and (c) preopen. 


Proor : (a) Suppose that f is quasi irresolute and semi-preopen. Let V € SO(Y). 
It follows from Lemma 2.5 that f~1(V) C s Int (f-1(s CI(V))). Moreover, by 
Theorem 5.5 we obtain f-! (V) C s Int (s Cl (f-! (V))). Therefore, it follows from 
Theorem 3.4 that fis almost irresolute. 


(b) Suppose that f is quasi irresolute and semi-open. Let V € SO(Y). By 
Lemmi 2.5, w2 hive f-!(V) Cs {nt (f-1(s Cl (V))). It is shown in Theorem 2 of 
Noiri!6 that fis semi-open if and only if f-1 (s Cl (B)) C Cl (f~ (B)) for every subset 
B of Y. Therefore, we obtain f-! (V) C s Int (Cl (f-! (V))) = Cl (Int (Cl (f-1 (V)))). 
By Theorem 3,4, fis almost irresolute. 


(c) Suppose that fis quasi irresolute and preopen. Let V € SO (Y). By Lemma 
2.5, we have f-1 (V) C s Int (f-1 (s C1(V))) C s Int (f-1 (CI (V))). Since f is preopen, 
by Theorem 11 of Rose2° we have 


f-1 (C1(V) = f-1 (Cl (Int (V)) C Cl (f=! (Int (V))) C C1 (Ff (VO). 


Therefore, we obtain f-! (V) C s Int (C1(f-1(V))) = Cl (Int (Cl (f-1 (V)))). By 
Theorem 3.4, f is almost irresolute. 


Remark 5.7 : The function f in Example 4.3 is almost irresolute, open and semi- 
preopen. However, it is neither quasi irresolute nor semi-continuous. 


6. STRONGLY SEMI-CLOSED GRAPHS 


For a function f: X > Y, the subset {(x, f(x)) |x € X} of the product space 
X x Yis called a graph of f and is denoted by G (f). Dube et al.9 defined and inves- 
tigated a strongly semi-closed graph. We shall improve on some results established by 
Dube et al... 


Definition 6.1\—The graph G (f) is said to be strongly semi-closed Dube® if for 
each, (x, ¥) & G (f), there exists U € SO (X,x) and V € SO (¥,y) such that [Uxs Cl 
(VV))]O G(S) = ¢. 

It follows immediately that G (f) is strongly semi-closed if and only if for each 
(x, y) & G (/f), there exist UE SO (X, x) and V € SO(Y, y) such that f (U) 
\ sCi(V) = ¢. The graph G (f) is said to be semi-closed (resp. closed) if it is semi- 
closed (resp. closed) in the product space X x Y. If G(f) is strongly semi-closed, 
then it is semi-closed. However, the converse is false by Example 2 of Dube et ai.®. 
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A function f: X¥ > Y is said to be pre-semi-open® if f (U) € SO (Y) for every 
U € SO (X). Every pre-semi-open function is both semi-open and semi-preopen. 
However, in Example 4.3, f 1 is semi-preopen but not semi-open and hence not pre- 
semi-open. Moreover, in Example 5.3, f is semi-open but not semi-preopen and 


hence not pre-semi-open. 


Theorem 6.2— Let f: X¥ ~ Y be either semi-preopen or semi-open. If G (f) is 
closed, then G (f) is strongly semi-closed. 


Proor : Since G (f) is closed, for each (x, y) & G(f), there exist open sets U 
and V containing x and y, respectively, such that f(U) N V = ¢. First, suppose that 
f is semi-preopen. Since V is open, we obtain s Int (s C1(f(U))) N sC1(V) = ¢. 
Since f is semi-preopen, f (U) € SPO (Y) and hence f(U) 1 sC1(V) = ¢ by Lemma 
3.2. This shows that G (/) is strongly semi-closed. Next, suppose that fis semi-open. 
Since f(U) € SO (Y), we have f(U) (1 s Cl (Vv) = ¢ and hence G(f) is strongly 
semi-closed. 


Corollary 6.39— If f: X > Y is pre-semi-open and G(/) is closed, then G (f) 
is strongly semi-closed. 


Noiri!8 showed that a semi-continuous function into a Hausdorff space has a 
semi-closed graph but it does not have a closed graph in general. 


Theorem 6.4— If f: X > Y is weakly quasi continuous and Y is Hausdorff, 
then G (/) is strongly semi-closed. 


Proor: For each (x, vy) & G(/f), there exist disjoint open sets V and W of Y 
containing y and f(x), respectively. Hence we have Cl (W) (4 Int (CI (V)) = ¢ and 
Cl(W) A s C1(V) = ¢. Since f is weakly quasi continuous, by Theorem 4.1 of Noiril® 
there exists U € SO(X, x) such that f(U) C Cl(W). Therefore, we obtain f(U) 
\ sCl(V) = ¢. This shows that G(/) is strongly semi-closed. 


Corollary 6.59— A function f: X + Y has a.strongly semi-closed graph if it 
has one property of the following : 


(a) f is semi-continuous and Y is Hausdorff; 
(b) f is weakly irresolute and Y is Urysohn; 


(c) f is 9-irresolute and Y is Urysohn. 


Proor: This follows immediately from Theorems 2.6 and 6.4 and Remark 4.2. 


A space X is said to be regular semi-T2 Maheshwari et al.1° (resp. semi-T2 
Maheshwari and Prasad!4 if for distinct points x, y of X, there exist disjoint regular 
semi-open (resp. semi-open) sets U and V such that x € Uand y € V. 


Lemma 6.6—A space X is regular semi-7» if and only if X is semi-7». 
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Proor : Every regular semi-7T2 space is obviously semi-7». Conversely, suppose 
that X is semi-T2. For distinct points x, y of X,there exist U € SO (X, x) and 
V € SO (X, y) such that U 1 V = ¢; hence UM s C1(V) = 6. By Lemma 2.2, we 
have s Cl (V) € SO (X) and hence s C1 (U) ™ s C1(V) = ¢. It follows from Lemmas 


2.1 and 2.2 that sCl(U) and s Cl (V) are regular semi-open. Therefore, X is regular 
semi-72. 


Theorem 6.7— If f: X > Yis quasi irresolute and Y is semi-T2, then G (/f) is 
strongly semi-closed. 


Proor : Let (x,y) & G(f). Since Y is semi-T2, there exist y € SO (Y, y) and 
W € SO (Y, f (x)) such thatV M7 W = ¢. By Lemma 2.2, sCl(V) 1 sCl(W) = ¢. 
Since f is quasi irresolute, there exists U € SO (X, x) such that f(U) C s Cl (W). 
Therefore, we obtain f(U) \ sCl(V) = ¢. This shows that G (/) is strongly semi- 
closed. 


A surjection f: X — Y is said to be semi-s-connected? if f~1 (V) is semi-clopen 
in X for every semi-clopen set V of Y. It follows from Lemma 2.5 that a surjection is 
set-s-connected if and only if it is quasi irresolute. Dube et al.9 defined a space X to 
be extremally-s-disconnected if s Cl(U) € SO (X) for every U € SO (X). However, 
by Lemma 2.2 the semi-closure of a semi-open set is always semi-open. 


Corollary 6.89— Let f: X — Y bea function and Y semi-TJ2. Then the follow- 
ing hold : 


(a) If f is irresolute, then G (/) is strongly semi-closed. 


(b) If fis a set-s-connected surjection and Y is extremally s-disconnected, then 
G (f) is strongly semi-closed. 


Proor: This is an immediate consequence of Remark 4.2 and Theorem 6.7. 


Corollary 6.92i— If f: X + Y is almost irresolute (in the sense of Thakur and 
Paik) and Y is regular semi-T2, then G (/) is semi-closed. 


Proor: By Theorem 2.6, “quasi irresolute’’ is equivalent to “almost irresolute” 
in the sense of Thakur and Paik2!. Therefore, this result is an immediate conse- 
quence of Lemma 6.6 and Theorem 6.7. 


Finally, we shall obtain a sufficient condition for a function to be quasi-irre- 
solute. A subset S of a space X is said to be s-closed relative to X, Di Maio and 
Noiri? if for every cover {V. | « € Vy} of S by semi-open sets of X, there exists a finite 
subset Yo of vy such thatS C U {sCl(Vu) | % © vol. If S = X, then the space X is 
said to be s-closed?. A space X is said to be extremally disconnected if the closure of 
every open set of X is open in X. 


Lemma 6.10— Let X be extremally disconnected and f: X — Y have a strongly 
semi-closed graph. If K is s-closed relative to Y, then f-1 (K) is semi-closed in X. 
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Proor: Let x & f-! (K). For each y € K, (x, y) & G (f) and hence there exist 
U(y) € SO(X, x) and V (y) € SO (Y, y)such that f(U(y)) Ns C1(V(y)) = $. 
Therefore, we have U (y) 1 f71 (s Cl (V (y)))= ¢ for each y € K. Since {Viy) | yEK} 
is a cover of K by semi-open sets of Y, there exists a finite number of points 1, 2 --., 
yn in K such that K C U {s C1(V (%)) | i = 1, 2, ..., m}. PutU = 1) {U (yi) | i=1, 2, 


_.., nN}. Since X is extremally disconnected, it follows from Theorem 2.9 of Jankovic! 
that U € SO (X, x). Moreover, we have U (\ f-! (K)=¢ and hence x & sCl ( f-1(K)). 
This implies that f-1 (K) is semi-closed in X. 


Theorem 6.11— Let X be extremally disconnected and Y s-closed. If f: X¥ > Y 
has a strongly semi-closed graph, then it is quasi irresolute. 


Proor: Let V € SO(Y). By Lemma 2.2, s Cl (V) is semi-clopen. It follows 
from Propositions 2.3 and 4.2 of Di Maio and Noiri’ that s Cl (V) is s-closed relative 
to Y. By Lemma 6.10, f~! (s Cl (V)) is semi-closed in X. Therefore, we obtain 
sCi(f-1(V)) C f-1(s Cl(V)). It follows from Lemma 2.5 that f is quasi irresolute. 


Corollary 6.12—Let X be extremally disconnected and Y s-closed semi-T2. A 
function f: X — Y is quasi irresolute if and only if G (/) is strongly semi-closed. 


ProoF: This is an immediate consequence of Theorem 6.7 and 6.11. 
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Three different sets of dual integral equations, involving Bessel function of 
order one, arising in some special axi-symmetric problems of elasticity 
theory and viscous flow theory are handled for their solutions, either in 
closed form or in terms of Fredholm integral equations of the second kind. 
The final results are in full agreement with the ones obtained through the 
Green’s function technique, used previously, for solving the mixed boundary 
value problems considered here. 


1. INTRODUCTION 

A number of workers have studied mixed boundary value problems, in the 
axisymmetric case, arising in Elasticity theory associated with punches and cracks!~3 
and in the theory of viscous flows induced by steady rotation or harmonic oscillation 
of a circular disc4-8. Of all possible methods of solution of these axi-symmetric pro- 
blems, the method of reduction to a set of dual integral equations!, appears to be the 
most natural and straightforward method of attack. But, as pointed out by Shail®” 
and demonstrated by Stallybrass®, certain dual integral equations are not amenable to 
their solution straightaway. It is for this difficulty that Stallybrass* and Shail*>°6 
have utilised an integral representation of the principal unknown potential by employ- 
ing a Green’s function technique devised specially for the purpose. 


We have shown in the present paper that by a suitable use of the Bessel equa- 
tion itself, it is possible to handle all the problems treated previously via the dual 
integral equations only and that the details of the Green’s function technique can be 

[oo] 
avoided, if we assume throughout the analysis that | p(s) Ji (sr) ds = noes | e-“ 
p(s) J1 (sr) ds, wherever such integrals occur. It must be emphasised that the Green’s 
function technique is better-suited to problems associated with more general axi- 
symmetric bodies than the circular disc for which the dual integral equations are the 
superior ones. 
2. THe Duat INTEGRAL EQUATIONS AND THEIR SOLUTIONS 
Problem | 
The dual integral equations 


Ts? A(s) Ju (sr) ds = — c(r),0<r <a ...(1) 
6 
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fey, (s) Ji (sr) ds = 0, r>a e8¢2) 
0 


arise (see Erguven®) in the study of a static penny-shaped crack problem in a homoge- 
neous isotropic elastic solid under torsion. Here o (r) represents the distribution of 
the shear-stress on the face of the crack and it is required that the displacement field 
given by the integral on the left of eqn. (1b) is zero at r = a, for the purpose of 
continuity. 


The method of solution of the equations (1) and (2) is as described below, and 
is different from the method described in Sneddon’s book!. 


We set 
Isa (s) Jy (sr) ds = f(r), O< r<a. ...(3) 


Then, using the well-known Hankel’s inversion formula to the eqns. (2) and (3), we 
obtain 


PGA j AS (A) Ja (As) dA. (4) 


The equations (1) and (4) finally give rise to the following integral equation for the 
unknown function / (A) : 


co a 
J s2 Jy (sr) ds f Af (A) Ja (SA) dA = —a(r), (0 <r < a) a be 
after interchanging the orders of integration, assuming that such an interchange is 
permissible here and even later on, for the other two problems treated in this paper. 
If we next use the idea that the Bessel function J (sr) satisfies the ordinary 


differential equation 


(+14 l 2 
as idee oe eae Ja (sr) = 0 ...(6) 


we observe that under certain special circumstances (as applicable to the crack pro- 
blem mentioned above), we can rewrite the integral equation (5) in the form 
ae 
. atl a ~1)u) =e), O<r<a) ..(7) 
with 


u(r) = H Ji (sr) ds | Af (A) Ji (As) dd. sehen 
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we then solve the ODE (7) by the standard method of variation of parameters and 
obtain the general solution in the form 


r 


+ 5 | o (t) dt — =| 126 (t) dt (9) 


oo 
0 2 
ee i pe 

r 


0 0 : ; , 
where c, andc, are arbitrary constants to be determined from the physical consi- 
derations of the problem. 


Of these two arbitrary constants, the constant a must be chosen to be zero in 


order that u is finite at r — 0, and the constant cy will be settled a little later on. 


Next, interchanging the orders of integration in eqn. (9), and using the result 
(see Shail5), 


co > min(r; A) 2d 

y y 
| Jy (sr) Jy (sA) ds -= | [(r2 — v2) (2 = v2)ft2 ...(10) 
0 0 


We obtain 


ae ee { fay dr 


lait | Q?— mie ~ ru(r), (0<r <4). .(11) 
v 


Using the Abel’s inversion formulae}, repeatedly to equation (1 1), we find that 


v 





C f(a) dr d 
r | iar = | an! (t))((v2 — ¢2) 1? dt, 


and, hence, 


+, fa une rm tae) ) ar 
py 


f(A)= 7a (a2—22)1/2% (a2 — ry) oe a — AQ2)1/2 ...(12) 


where 





ee aah) 
b (2 = myn 
Equations (12) and (13), along with eqn. (9) completely solve the ak equation 
(5), if the constant cy is determined. In order to determine the constant c cy , we need 


to use the result (12) and the observation that f (4) = . 0, arising out of the physical 
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requirement involving the continuity of f(r) atr =a, as argued earlier. We then 
obtain the equation that 


a da 
f at (tu (t)) dt 5 (14) 
6 (a2 — 1212 “y 


and this, together with eqn. (9) serves as the determining equation for the constant 


0 


Cc 


As a particular case, if we take o(r) = cr, wherec is a known constant, as 
considered by Erguven!9, we easily find that 





u(r) = cy + a r2 CLS 
and eqn. (14) gives that 

oe za ...(16) 
so that eqn. (12) decides that 

[A= “s A (a2 — A2)1/2 StF 


which agrees with the result of Erguven?!®. 


Problem 2 


The following dual integral equations arise in the study of the dynamic Reissner- 
Sagoci problem considered by Shail® : 


- a ge 
fe@A Gedi= “<p <i) (18) 
J (E+ p?/B2)1” C (&) Ja EP) dé = 0, (P > 1), (19) 


In order to derive the above dual equations, we have used the representation of the 
solution of the mixed boundary value problem Of Shail®, in the form : 


0 (P, 2) =f CE) e@ +8 1x P20 
12) =P C@e Hebe) dt, {> ot ...(20) 


which must satisfy the p.d.e. and the boundary conditions as given by 


02 0 1 00 7 02 a 27 
oS Oh 5 po 
ép2 * P ep p2 Bars en are: 0,z>0 oee(21) 
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and 
_ ap 
aah, | Samey 0 S =< } 
p= 5 .0<P<)) 
Fs on z = 0. a 64 
ee = 0,(P > 1) 


Shail® has used, like Stallybrass9, a special Green’s function technique to solve the 
mixed problems (21) and (22), as no direct method of attack exists to solve the dual 
eqns. (18) and (19). We present an approach, which is similar to the one employed 
for Problem | above, utilizing the Bessel equation (6), and show that for large values 
of p, the solution of the dual equations (18) and (19) can be determined exactly in the 
same manner as described by Shail®. 


We set 
c@-(5+ 2 yp" D® ...(23) 


and rewrite eqns. (18) and (19) in the form 





ie + a D (&) Ja (EP)dE= = 0 <p <1) (23) 
0 
and 
[(S+e ) D (&) Ji (EP) dé = 0(P > 1). (25) 
0 
Equation (25) can be recast in the form 
Ge Li aH aa =~ fe) Pas (Ep) dE = 0, (P > 1) 
...(26) 


after using the Bessel equation (6), and this, on integration gives that 
f D(@) Jy (Ep) d& = Co Ki (pP/P), (P > 1) -Q7) 
0 


after neglecting the part of the solution involving T, (pe/B), where I; (x) and Ky (x) 
are the modified Bessel functions of the first order and Co is an arbitrary constant to 
be determined by using order physico-mathematical consideration. 


We thus find that the dual equations (24) and (25) can be recast into the form 
(24) and (27) respectively. 
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If we then set 
°(p3yB2 + ERE D (E) Ja (Ep) dé = 8 (P), (P > 1) _..(28) 
0 
and use the Hankel inversion formula to eqns. (24) and (28), we obtain that 


(pip? + PDO = he +k | Ag (A) Ja (EA) da (29) 
nk 


(See Gradshteyn and Ryzhik!!, p. 683) and that eqn. (27) finally gives rise to the 
following integral equation for the unknown function g (P) : 


f p a { Jo (E) Ja (EP) 
1 0 
(Pp > 1) ...(30) 
where 
as aa . 
MP, = | eeagayre 11 PY i rd dé an 
0 


The integral eqation (30) is similar to the one obtained by Shail®, and can be at- 
tacked for solution, for large p, by a repeated use of the Abel’s inversion formulae as 


described by Shail®, if the following asymptotic result for the kernel M (P, A) is made 
use of : 


co 
ee +r e~ 2aw dw . 
M (P, A) = (PA)L2 ea(pta) | rm (w—p)1I2 (w—a)il2 +O (q 2), 
max(p 
; ...(32) 
with gq = p/p. 


We do not proceed any further with this problem here, as the other details are 
going to be repetitions of Shail’s work®. 


Problem 3 
The following dual integral equations arise, as shown by Goodrich’, in the study 


of a viscous flow problem induced by a rotating circular disc, kept on the surface of a 


bulk fluid of viscosity », which is otherwise contaminated by an adsorbed fluid film of 
different viscosity ». 


r f(y) Ji (yr) dy = wr, (0 <r Sa) ...(33) 
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F (uy + ny) f(y) Ja (yr) dy = 0, (r > a) ...(34) 


where “ is the constant angular velocity of the disc of radius a, rotating around its 
axis. 


Goodrich‘ has devised a special method of solving the above dual equations (33) 
and (34) in the three different circumstances, as given by the cases (i) » = 0, (ii) »=0 
and (iii) » #0,» 40. But, as pointed out by Shail5, Goodrich’s solutions are not 
the correct ones, since they involve certain divergent integrals. 


{t is because of this major difficulty that Shail5 has attacked the physical pro- 
blem of Goodrich? and several other generalization of it®-8, by a method utilizing the 
Green’s function technique. 


We have shown below that the dual integral equations (33) and (34) can also be 
attacked for their solution in the three cases (i), (ii) and (iii) as considered by Good- 
rich4, in a straightforward manner as described for the Problems | and 2 above, and 
we thus infer that the use of the Green’s function technique can be avoided here 
also. 


Case (i): » = O— In this special case, the dual equations (33) and (34) take up 
the forms : 


F£0) A On) dy = or (0S <a) (35) 
Py? £0) Ji (yr) dy = O(r > a) ‘ ...(36) 


and, the second eqn. (36), can be recast, by using the Bessel equation (I/), in the 
form 


if £0) J, (yr) dy = Colr (r > a) ...(37) 


where Cp is an arbitrary constant to be determined. 


A straightforward use of the Hankel inversion formula to eqn. (35) and (37) 
gives 


St (”) = Co Jo (ay) + wa? Jo (ay). (38) 
(see Gradshteyn and Ryzhik!!, p. 683). 


We ultimately find that the constant Co appearing in (38) must be chosen to be 
wa” in order to make the integral in (36) convergent, in the sense mentioned in the 
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introduction, and this agrees with the observation of Goodrich, even though f (¥) is 
different. 


Case (ii) : 1 = O— In this case, the dual equations to be be solved are the ones 
as given by eqn. (35) and the new equation 


f »f 0) Aa On) dy = 0, (r > a). -..39) 
Assuming that 
Por 0) On) dy = 8 (1), (0 <r <a) (40) 
and using Hankel’s inversion formula we find that 
oo 
f(y) = J g (A) A Jy (Ay) ay. ...(41) 


Then using (41) in eqn. (35) we ultimately derive that 


r a 

2 | v2 dv g (A) da 

rr (r2 as y2)1/2 (A2 —_ y2)1/2 
0 v 


= or, (0 <r « a) ...(42) 


obtained after utilizing the formula (10). 


A repeated Abel inversion procedure, like the one adopted in the previous pro- 
blems, ultimately gives 


m g(r) = 4ur (a2 — r2)-12, (0 < r Ea) ...(43) 
and the solution of the dual equations can be completed by using (41). 


Case (iii): » # 0, 1  O— In the most general case of the dual equations (33) 
and (34), when neither » nor 4 is zero, we obtain a Fredholm integral equation of the 
second kind, which is similar to the one obtained by Shail5, by means of a procedure 
as described below. 


We first rewrite eqns. (33) and (34) in the form 
oo 
Lf) Ji (yr) dy = or O <r <a) ...(44) 
and 


r y (1 + Ao ay) f(y) nr) dy = 0 (r > a) (45) 
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where Ap = »/ua. 


Then, setting as in the previous problems, 
io 2] 
J y(l + Asay) f(y) Ii (yr) dy=g(r) O<r<a) ...(46) 


and using Hankel’s inversion formula to eqns. (46) and (45), we find that 


a 


1 
(y= Trica | A g (A) Ja (VA) A. ...(47) 


We next use the relation (47) in (44) and interchange the orders of integration to 
obtain the following integral equation for the function g (r) : 


a co 


| Ag(A) da \(es) Jy (yr) Ja (y)) dy = or, (0 Sr <a). ...(48) 
0 0 


Equation (48) is an integral equation of the first kind and it can be converted 
into an integral equation of the second kind, by observing first that 


1 | 


iia Cee ESE gor: 


and then using the result (10) along with an Abel’s inversion procedure, of a type 
similar to what is known as Williams’s method!2. 


We find that the second kind Fredholm equation is obtained finally in the form 


8s A Kl f L* (t, v) g* (t) dt = 4wy, (0 < v < a) ...(50) 
0 
where 
r (A) da 
g* (v) = vl ce cok 7 W 
and 
3 oe We 
Ae ee A OY : : Glee 
L* (t,v) z | —pe sin (ty) sin (vy) dy (52) 
0 


Equation (50) can be easily identified to be similar to the one obtained by Shail®, for 
the problem of Goodrich’. 
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A slightly more general dual integral equations, than the ones given by (33) and 
(34), which are 


if ai rs aa aeay ai Jro Ji (yP) dy = op 0 < P <a) ___...(53) 
and 
2 ro) y Ji (vP) dy = 0, (P > a) ..(54) 


where h > 0, can also be reduced to a Fredholm integral equation of the second kind 
by a method similar to the one described above. 


The equations (53) and (£4) arise in the rotating disc problem of Shail®, when 
the disc is kept at a distance a below the contaminated surface considered by Good- 
rich4, so that the particular case h = 0 of (53) and (54) correspond to eqns. (33) and 
(34). 


3. CONCLUSION 


The present paper has dealt with the dual-integral-equations-formulation of some 
of the well-studied axi-symmetric mixed boundary value problems of potential theory 
in a manner different from the ones used in the literature, but is very straightforward 
otherwise. The principal aim has been to show that for the three problems discussed 
here, or its generalizations, it is not necessary to employ any other complicated 
technique, such as the Green’s function technique, used by previous workers even 
though the merit of the latter technique is unquestionably of a higher level, if one has 
to handle axisymmetric bodies other than Circular discs as has been the case with the 
above three problems. 
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BANACH SPACE VALUED DISTRIBUTIONAL MELLIN 
TRANSFORM AND FORM INVARIANT LINEAR FILTERING 
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Department of Mathematics, Regional College of Engineering, Nagpur 440011 


(Received 24 June 1987; after revision 17 March 1988) 


From invariant filters are those shift variant filters such that a linear 
scaling of their input gives rise io a linear scaling of their output. In this 
work we develop a theory of Mellin transform and Mellin type convolu- 
tion for Banach space valued distributions. Application of this theory to 
relate the input, output singals of the form invariant systems is given. 


1, INTRODUCTION 


Any physical system, be it electrical or optical, can be characterized by an input 
signal, output signal and system function. Time invariance (shift invariance) is a 
property possessed by many physical systems. Zemanian!” introduced the concept of 
Banach space valued distributions. He presented a theory for the convolution and 
Laplace transformation of Banach space valued distributions and used these concepts 
for convolution representation of continuous linear time invariant systems. 


The form invariant filters are those filters, such that a linear scaling of their input 
gives rise to a linear scaling (generally through a different scaling factor) of their out- 
put. Braccini and Gambardella®’ defined the form invariance property and obtained 
the most general class of linear shift variant systems satisfying it. They related the in- 
put and output signals of the form invariant systems by using Mellin transform. 
Further they discussed applications of form invariant filters in optical pattern reco- 
gnition, image restoration, image reconstruction, radar signal processing etc. Braccini 
and Gambardella4 considered the most general input output relationship for a linear 
one dimensional system as 


g(x) =_f f(s) w(t, x) de 


where f(t) is the input signal, g (x) is the output signal, and w (#, x) is the kernel of 
the integral transform. The integral representation has no sense when the input signal 
f (t) is a singularity signal. For example singularity signal like 5 cannot be treated as 
Lebesgue integrable function (LIF). If, we treat 5 as LIF we get senseless results (see 
Gelfand and Shilov®, p. 4, Zemanian®, p. 10). The theory of distributions (generalized 
functions) provides a natural language for such signals (see Zemanian’, Wohlers and 
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Beltrami’). By using distribution theory we can generalize and correct results of 
of Braccini and Gambardella‘. 


Recently the author9-11 extended some integral transforms to Banach space 
valued distributions and discussed their applications to system theory. 


In this paper we develop a theory of Mellin transform and Mellin type convolu- 
tion for Banach space valued distributions. We use this theory to relate the input, 
output signals of the form invariant systems. Justification for these studies includes the 
desire to admit a larger class of systems and input output pairs. When the input signal 
f(t) is a LIF our results of section 5 reduce to results of Braccini and Gambardella?. 


The plan of the paper is as follows. In section 3, we develop a theory for Mellin 
transform of Banach space valued distributions. In section 4 we present a theory for 
Mellin type convolution of Banach space valued distributions. Finally in section 5, 
we discuss Banach space valued distributional form invariant linear system and use 
Mellin transform to relate the input output signals of the form invariant systems. 


2. NoTATIONS AND TERMINOLOGY 


The notations and terminology of this work wil] follow that of Tiwari91011 and 
Zemanian!, Rand C denote, respectively, the real line and complex plane R, de- 
notes the positive half-line {f € R:0 << t < cco}. A and B denote complex Banach 
spaces. E (A) denotes the linear space of all smooth functions ¢ from R into A. E, 
(A) is linear space of all smooth functions from R, into A. Let K bea compact sub- 
set of R. Dx (A) denotes the linear space of all smooth, A valued functions ¢@ such 
that supp (¢) C K. The space Dx (A) is assigned the topology generated by the collec- 


tion {rk} °°, of seminorms, where 


rk (¢) = sup |[d* (t)ll4. 
tEK 


| . ||4 denotes the norm in A. Let {Ki} be a sequence compact subsets in R such 


oo ~ 
that Ky C°Rs.C fs Y Kj = R. Thespace D(A) is defined as the inductive limit 
space generated by Dx, (A). That is 


D(A) = ind Dg (A). 
aa 


Similarly we have 
D,(A) = ind Dr, (A) 


ie 


here Ki, Ke..., are compact subsets of Ry, and R, = 
3 


“C8 
a 
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For A = C, we write E (A) = E, E, (A) = E,, D(A) = Dand D, (A) = Dy. 
When U and V are topological spaces, [U;V] denotes the space of all continuous 
linear mappings of UintoV. The symbol </, ¢> denotes the element of V assigned 
tod € Uby f E [U; V]. The notation [1 denotes the end of a proof. 


3. MBLLIN TRANSFORMATION 
Following Tiwari!2 we define the space Ma,v,. (A). 


3.1. The Space 
May,b,« (A) 
For « > 0, we define 


May,b,« (A) = {:¢ € E, (A), 


ia,,k ($) = sup ||Aao (t) #72 DY, ¢ (2) 
Ry 


= 1CEL* ke, k == 0, 1,....}. 
The constants and Cx depend on the function ¢ and 
en fey.) 
Rare CS = lies $< .00: 
For k = 0, weset k** = 1. The topology of the space Ma,v,« (A) is generated by 
the family of seminorms '{ia,v,t};° - The space Ma,» (A) is defined as the inductive 
limit space generated by Ma,p,« (A). That is 


Ma,» (A) = ind May, (A). 


a—->oo 
Following Zemanian! (p. 102) we define the space M (w, z; A) as 
M (w, z; A) = mes bd, (A) 
where 
an >w,, bn > 2_, Ww, z € [— ©9, ©]. 


Let B be a complex Banach space. Any f € [M (™, 2; A), Blis a Banach space — 
valued distribution. When 4 = B= G, f becomes a scalar valued distribution. 


Let S be a collection of bounded subsets of Ma,» (A). The topology of uniform 
convergence on [Ma,» (A), Bis that generated by the collection of seminorms 


{s where, 
Koes 


ca(f) = sup || <f,¢ > Ile. 
$EQ 
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Se [Ma,» (A), B}, QE S. { 


The weak topology or simple topology for [Ma,» (A), B] is generated by the colle- 
ction seminorms {Pg}, where ¢ traverses Ma,» (A) and 


Pe(f) =< fh¢> Ib. 
We now define Mellin transform of A-valued and [A, B] valued distributions (see 


Zemanian!, p. 115). 


3.2. MBLLIN TRANSFORM OF A-VALUED DISTRIBUTIONS 


Let f € [D,; A]: We way that f is Mellin transformable if there exists two 
members oj, cg € [— oe, cc] such that «1 < oo, f € [M (a4, 62): A], andin addition 
f & (Mv, z); A] if either w <0, or z > og. With Of = {s: 01 < Re(s) < 939}, 
18-1 € M (01, 52) we define Mellin transform Mf of f as 


F (s) =(M f)(s) = < f(t), #-1>, 5 € ay. staet) 
It can be easily proved that F (s) is an A-valued analytic function on Qy. 


The space [M (w, z; A); B] can be identified with the space [M (w, z); (A; B)| 
through the equation : 


<fy ¢>a= <y, ga> 
(see Zemanian!, p. 105) where f, € [M (w, z); [A; B]], 
y € [M (wv, z; A); B],d © M(w,z) andae A. 


Because of the above identification we use the same symbol to denote both fy 
and y, and define Mellin transform of [A; B) valued distribution as : 


We say that y € [D,;: A] is Mellin transformable if there exists 1, %2 € [-°, 
eo] such that 11 < y2,» € [M (m1, y2; A); B] and y & [M (¥,z; A); B]if either w < yy 
or Z > 72. 


Using the above identification y € [M (1, 2); [A; B]] also. Hence we now de- 
fine the Mellin transform Y of y as 


Y (s) = <y(t); #1>,s EQ, met 57) 
where %y = <y(1); -1>,5 € Q, (3.2) 


where 0, = {s:y1 < Re(s) < yp} is called strip of definition for the Mellin transform 
of y. 


Theorem 3.1 (Analyticity Theorem)—If My = Y(s) fors € Q,, then Y (s) is an 
[A; B) valued analytic function and for each non negative integer k 


YF (s) = <y(1), DY t-1> 
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PRoor : From definition the result is true for k = 0. With fixed s and As + 0. 
consider 


<y (t), Yas (t)> 


Y* (s + As) — Y* (s) bone 
Aiea Sa — <y (t), D, sot SP ae 


It is not difficult to show that vas (rt) converges to zero in M (n1, 2) as As > 0. 


4. MELLIN TyPpE CONVOLUTION 


Theorem 1.1—The scaling operator Sa: ¢(t) > ¢(at),a > Oisa topological 
automorphism on the space Ma,»,. (A). 


Proor : Sq is clearly well defined and linear. For continuity we observe that 


sup l| Aa (t) t**2 Dr ¢ (a. 2)I| 


I 


ae ll Aa, 6 (=) (> yen at D ¢(T)|, T = at 


C sup ||Aa,» (T) T#t1 D* ¢ (T)||, Cis a constant 
R 


+ 


=C la,b,k (¢). 
: -1. aia ft 
The inverse mapping S_ is defined by S>:¢(t) ~¢ ( - ) . 
Theorem 4.2—If ¢ € M (w, z), then, for each fixed ¢, P(t x) E€ M(W, z)asa 
function of x, where x > 0. 


Proor: Proof follows immediately from Theorem 4.1. 


Theorem 4.3—If f is a member of [M (w, z); A] and ¢ € M(w, z), then d > b 
is a continuous linear mapping of M (w, z) into M (w, z; A), where 


b(t) = < f (x), ¢ (t x)>. (4.1) 


Proor :We first prove that  (t) is an A-valued smooth function on 0 < f < oo, 
i. e. we want to prove 


UE (t) = < f(x), D) o(tx) >. (4.2) 


The above result is true for k =~ 0 by definition. 


Fork = 1, ¢t fixed and Ar # 0 consider 


Het ana O _ <f (x), Dig (tx)> = <S (x), Oat (x)> -..(4.3) 
t 
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where 
doe (x) = ge [6 M(t + At) x} — 4 (Ew) — Did (>). 


We now show that 6az (x) converges in M (w, z) to zero as At — 0. Using Taylor’s 
formula and treating At as independent variable, it can be easily proved that 


| Aa,o (x) xPtl 1 (x)} | + 0 as At — 0. 
Thus result is established for k = 1. Now by induction the proof can be easily com- 
pleted. 
Next to prove ¥ (t) € M (w, z; A), using the boundedness property of distribu- 


tions we note that 


Aan (t) t#*2 DE o (2) 
= [Aa (t) #41 DF < f(x), 6 (tx) > | 


<M sup | Aan (t) tt Aaso (x) xP*1 D? Di $ (tx) | 
R 


+r 


< M sum | Aa,o (t x) (t x)Ptk+1 DP*k g (t x) | 
R 


+ 
S M ia,v,p..k ($). ...(4.4) 
Thus } (t) € M(w, z; A). Continuity of the mapping ¢ > w also follows from (4.4). 
Theorem 4.4—The mapping f — defined by 
¥(t) = <f (x), ¢(t x)> 
is a linear mapping that is uniformly continuous with respect to S sets in M (w, z). 


ProoF : Let Q be any S set in M (w, z). We choose a and b such that w < a<bh 
<z. Then Q is a bounded set of Ma,». Now U* (t) = <f(x), 6* (tx)> 
so that 


Aayo (x) x”*1 D™ [ayo (t) +1 DF ¢ (2 x)] 
== Aa,p (t x) (t x)Ftmtl Dk+m ¢ (¢ x), -+.(4.5) 
From (4.5), for k fixed, as ¢ traverses Q and f traverses Ry, Aa, (t) t*+1 D: ¢ (t x) 


as a function of x traverses a bounded set in @ in Ma,. This means that Ois aS set 
in M (w, z), and we have 
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sup sup ||Aa,o (t) t**1 UL (£)I] 4 


gEQ tERy 
= sup sup || Aa,v (t) t**1 < f(x), pk (t x)>Kl4 
¢E€Q tER, 
= sup |< f, @> lla. ...(4.6) 
GEO 


Now if A is a neighbourhood of zero in M (w, z; A), (A) (\ Ma,» (A) C (A) isa 
neighbourhood of zero in Ma,» (A). From (4.6), there exists a neighbourhood = of 
zero in [M (w, z); A] such that | € (A) for all dé € QO andallf € &. 


We now define Mellin type convolution. 

Mellin Type Convolution—Let a and b be two real numbers with a < 5. Mellin 
type convolution is an operation that assigns to each arbitrary choice of the pair 
y € [M (y, z; A): B] and f € [M (w, z); A] the product y v f defined by 

<yvf,o> = <y(t), <f (x), 6 (tx)>> 
¢ € Mw, z) (4.7) 

Theorem 4.5—Assuming the validity of Theorems 4.2, 4.3 and 4.4. The operator 

yyv:f ~ yy/fis a continuous linear mapping of [M (w, z), A] into [M(w, z); B]. 


Proor : Linearity is clear from (4.7). To prove continuity we have to prove 
that any neighbourhood A of zero in [M (w, z); B] contain a neighbourhood & of zero 
in [M (w, z); A]. 

Let ® be an arbitrary s — set in M(w, z). We have 


oo(yvf)= sup || <yvf, o>llB 
E0 
= sup || <y, ¥> lle. 
¢E® 


Because y € [M (w, z; A); B], to each neighbourhood I of zeroin B there exists a 
neighbourhood © of zero in M (w, z; A) such that y maps Q into [. Now the proof 
can be completed by using the fact that y v fis a composite mapping ¢ > / > <y, U>. 


We now prove in the following theorem that the operator y v commutes with the 
scaling operator Sz. 


Theorem 4.6—If a > 0, y € [D, (A) Bl and f € [E,, A], then 
Sa(yvf) = y (Saf). 


. Proor : By Theorem 4.1 Sa: ¢ (t) > ¢ (a4) is an automorphism on D, and 
hence 


<Sa (vf (> = <OS) = (+= )> 


(equation continued on p. 500) 
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=<y()< sf), + 3 (<)>. (4.8) 


a 


Also 
= v (Sa f ), C) (t x)> 


= <y(t) <Saf (x),  (t x)> 


= <y(t) <f (x), ~ #( Fa )>. (4.9) 


a 


From (4.8) and (4.9) our theorem is proved. 


Theorem 4.7—If f © [Ma,»o (A), B] andw € D, (A). Then, w > gisa conti- 
nuous linear mapping of D, (A) into E, (B), where 


g(x) = <s0, 2 w (*)>. 


ProoF: Itis easy to prove that g is smooth and the mapping is linear. To prove 
continuity we observe that 


lig* (lle = | < f(t), D* | ~w(+ )> 


~ 


< max sup ||D. Di - w (=n 
O<l<r R, : 


From above continuity easily follows. 
We call 
] 
(fri = <f(), > w( x )>. 


The Mellin type regularization of f by w. 


Theorem 4.8—If y € [D, (A); B] and My = Y(s) fors € Qy, if f € [D4; A] 
and Mf = F(s)fors € Qs, and if Qy 1) Qs is not empty then y v f exists in the sense 
of Mellin type. Convolution in [M (w, z), B] where the interval (w, z) is the intersec- 
tion of Q, M OF with the real axis. Moreover, 


M (yv f) = Y(s) F(s). 


Proor : It is already proved in Theorem 4.5 that yvf € [M (wv, z); B). Further 
because 1*-1 € M (w, z) for each fixed s with w < Re (s) <2, 


M(yvf) = <y(t), f(x), (tx)1> 
= <y (t), #-1> <f (x) xt-1> 
= Y (s) F(s). 
Note that Y (s) F(s) is B-valued function for any fixed s € Qy 1 Qy. 
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5. FORM INVARIANT LINEAR FILTERING 


Generalizing the result of Braccini and Gambardella* (p. 1613) we write the in- 
put output relationship for a linear one dimensional system as 


g(x) = <f(t), w(t, x)> (5.1) 


where f (t) € [E (A), B], w(t, x) © E (A) and g (x) is the output signal. When A 
= B = C the Banach space valued distribution f becomes scalar valued distribution. 
Further if f(t) is a regular scalar valued distribution generated by locally integrable 
function f, we can write (5.1) as 


co 
g (x)= bh f (t) w (t, x) dt. (52) 
Further generalizing the definition of Braccini and Gambardella4 we define the Banach 


space valued distributional form invariance property as below : 


Let ga (x) be the output of the system (5.1) when the input f (t) is replaced by 
f (at), a being a positive real number. We say that the system is Banach space valued 
distributional form invariant if and only if 


Za (x) = « g (8x) 33(535) 
where « and 8 are real functions of a. From (5.1) 
ga(x) = <f(at), wt, x)> 


and 


ag(Bx) =a <f(t), w(t, 6B x)>. 
From (5.3) we get 
<f (at), w(t,x)> =« <f (t), w(t, Bx)>. (5.4) 


Using the following property of distributions, namely 


<fat,¢@> = <f, } 6(Z)> 
we get from (5.4) 
- <f(t), w( 4, x)> =a <f(t), w(t, Bx)>. 


The above equation is true if and and only if 


ae w( oy, x) = «w(t, Bx). .. (5.5) 


a a 


General solution of this equation is (see Braccini and Gambardella4, p. 1613). 


w (t, x) = x°® Wo (t/x?) (5.6) 


B=alle, a= a(s-e)'e, x #0 
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where 5 and c are real numbers. In (5.6) choice of o affects the scale factor B of the 
output signal « g (8 x), whereas the choice of 5 (for any given value of c) affects the 
amplitude factor « of the output « g (6 x). 


We now discuss some particular cases of (5.6) and relate the input output signals 
of the form invariant systems using Banach space valued distributional Mellin trans- 
form. 


Case !—Taking 5 =o = 1, B =a, « = 1 in (5.6), we get 
w (t, x) = x-1 wo (t—x), x € 0. Sk be BP i 
After some simple manipulation (5.7) can be written equivalently as 
w (t, x) = t71 wo (x/t), t + 0. s<( Do) 
Note that taking 8 = a implies the same scaling factor in both the input f (at) and the 
output « g (8 x) of system (5.1). Further « = | means no amplitude gain changes are 
undergone by the output « g (Bx), when the input is linearly scaled. 
Now using (5.8) we write input output relationship from (5.1) as 
g(x) = < f(t), t-1 wo (x/t)> = (fv wo). ...(5.9) 


Note that here f is a Banach space valued distribution and not an Ordinary integra- 
ble function as in Braccini and Gambardella. 


Taking Mellin transform of both sides of (5.9), we get 


G(s) = F(s). Wo (s). .--(5.10) 


(5.10) is similar to relationship relating to the Laplace transforms of the output and 
the input Of a shift invariant (time invariant) systems (see Zemanian!). The result is 
in agreement with result of Braccini and Gambardella4 (p. 1618) when f(t) is a 
Lebesgue integrable function. 


To discuss next particular case of (5.6) we need 


Theorem 5.1—The operator P : Ma-r,0_r1,0 (A) > Ma,v,. (A) defined by P (¢) 
= 1 ¢ is an isomorphism from the space Mayr,o_r,« (A) onto Ma,»,. 


real number. (A) whereris a 


Proor : Observe that 


k 
af lAa,v (t) t*+1 D, (t" $)ll 


k 
= sup |lAa,» (t) t#+1 SZ Cp pr-v ja tee 
R, on Dp t $|| 
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where Cp is some constant. 


k kp 
= sup || rai be Aa-r, v-rt**1-P Dll 
p= t 


S 
r 


Cp la-r,b-1,k-p (¢). 


i 
LMer 


This proves the continuity of the map P. Linearity of the map is easily seen. The in- 
verse mapping p~1: Ma,v,. (A) > Ma-_r,v-1,2 (A) is defined by p-! (¢) = 1-7 ¢. The 
linearity and continuity of p~! can also be proved similarly. Thus P is an isomorphism. 


Case Ji—Taking w (t, x) = t" x" wo (x/t) where » andy are arbitrary real 
numbers. 


21x) =< 7 (1), w(t, x) > ad 


I 


g (x) = < f(t),.t-1 11 xP wo (x/t) > 


[ f(t) v t%*1 x* wo (x/t)). 


Using Theorem 5.1 and some simple properties of Mellin transform (see Snedonn}’, 
p. 270) we have 


G(s) = F(s+yt+ut 1) wo(s + p). 


This result is generalization of the result obtained by Braccini and Gambardella* 
(p. 1618). 


CONCLUDING REMARKS 


In this paper we have applied Banach space valued distributional Mellin trans- 
form to 1-dimensional linear filters. In future we plan to extend these results to 2- 
dimensional linear filters. For this we will develop two dimensional Mellin transform 
and Laplace-Mellin (or Fourier-Mellin) transform. These transforms will be used to 
relate the input output signals. We will also prove that form invariance implies the 
operator has a convolution (Mellin type) representation. 
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A CYLINDRICAL WAVE-MAKER IN LIQUID OF FINITE 
DEPTH WITH AN INERTIAL SURFACE 
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This paper is concerned with the generation of waves ina liquid of uniform 
finite depth with an inertial surface composed of uniformly distributed non- 
interacting floating particles, due to forced symmetric motion prescribed on 
the surface of an immersed circular vertical wave-maker. The techniques of 
Laplace transform in time and a suitable Weber transform in the radial 
co-ordinate are used to solve the problem. It is shown that if the inertial 
surface is too heavy, time-harmonic disturbance due to the wave-maker has 
only local influence on the liquid. 


1. INTRODUCTION 


The classical wave-maker problem for the case of deep water with a free surface 
was solved long back by Havelock? wherein the wave-maker is either in the form of a 
vertical plane or a cylinder with circular cross-section. Later Rhodes-Robinson? 
extended them to include the effect of surface tension at the free surface. Recently 
there has been a considerable interest in different problems concerning generation of 
waves in a liquid with an inertial surface composed of athin but uniform distribution 
of disconnected particles (e.g. broken ice, floating mat on water). Rhodes-Robinson$, 
Mandal and Kundu2>, Mandal? considered problems involving generation of waves 
at an inertial surface due to different types of sources with time dependent strengths 
submerged in a fluid of both infinite and uniform finite depths. Rhodes-Robinson® 
also pointed out briefly the method of solving the plane-vertical wave-maker problem 
in a liquid with an inertial surface by a suitable use of Green’s integral theorem in the 
liquid region after taking Laplace transform in time. However the circular cylindrical 
wave-maker problem needs attention as it can, but not easily be solved by this method. 
Here we use a suitable Weber transform! in the radial coordinate after employing 
Laplace transform in time to solve the problem. The important time-harmonic case is 
considered and the inertial surface depression is calculated. It is observed that the 
time-harmonic wave-maker affects the inertial surface only locally if it is too heavy. 


2. STATEMENT AND FORMULATION OF THE PROBLEM 


We consider the motion under gravity in an ideal liquid of volume density P 


covered by an inertial surface composed of uniformly distributed floating particles of 
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area density Pe. ¢ = 0 corresponds to a liquid with a free surface. On an immersed 
vertical circular cylindrical wave-maker, the normal fluid velocity is supposed to be 
prescribed which is both time and depth dependent. We choose a cylindrical coordi- 
nate system (r, @, y) in which the y-axis is taken as the axis of the cylinder with radius 
a so that r = a is the wave-maker0 <<» <Ar>a isthe fluid region and y = 0, 
r >a is the position of the inertial surface at rest. The wave-maker starts operating 
from time t = 0 with outward normal velocity U(y, t) on its boundary r = a. We 
consider only the axisymetric case in which the resulting motion in the liquid is 
independent of 6. Since the motion starts from rest it is irrotational and can be 
described by a velocity potential 9 (r, y, t) satisfying the Laplace’s equation 


| 
Ger Or Bey 80,6 > By Ory ae em 35 8 
The condition at the inertia) surface y = 0 relating the potential function and the 


inertial surface depression ¢, within the frame-work of linearised theory, consists of 
the kinematic condition 


ey = t ony =0 i gies 
and the dynamic condition 

ge =ett+e lu wosl2ea? 
where g is the gravity. Elimination of produces the inertial surface condition 

On —go,=0ony=0 .. (2.4) 
where 

® = 9 — 29, ...(2.5) 
The condition at the wave-maker is 

g = U(y,t) onr=a .+.(2.6) 
and the condition at the bottom is 

gy =~O0ony=A, ea ( 2am) 
There are also initial conditions at the inertial surface given by 

Lae =@=0 on y=0 at t=0, Vaal Qiao) 


Let a bar above a function denote its Laplace transform in time. Then @(r, y; p) 
satisfies the boundary value problem 


= ] 
Grr + > Reon yy = Cnr > eee 


P?o — (g +€p2)3,=0 on y=0,r>a 
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or = U(y,p) on r=a 


SG =90 on y= A. ee eh) 


3. SOLUTION BY WEBER-TRANSFORM METHOD 


We use the following form of Weber-transform of a function f(r) defined in 
(a, cc) by 


g(t)= | rA(r,®) f(r) dr & > 0) 8.1) 


where 


A (r, &) = Ja (@&) Yo (r &) — Jo(r &) %1 (a 8) Asie) 


Jn, ¥n (n = 0, 1) are the Bessel functions of the first and second kinds respectively. 
(cf. Davies!, p. 252). 


The inverse transform formula is 


fi) = | Be ANG) og (5) at. 3.3) 
2 J? (a) + Yi (a8) 


It may be noted that 


oo 
] 2 
| ( fir + ca fr )r A (r, e) d= —- as fr (a) ones E2 g (&). ...(3.4) 
Let & (&, y; p) denote the Weber transform of & (r, ¥; p) as defined by (3.1). Then 
in view of (3.4), v (r, &) satisfies 
es } 
bw 89 = = T0;p), 0<¥<h 
| ealae) 
ped —(g+ep)ty=0 aty=0 | 
by = 0 at wi = h. J 
The solution of (3.5) is 
..3.6) 


h 
be xp)= — 2 | 60,2) O(%p) ds 
0 


where G (y, «) is the associated Green’s function given by (cf. Mikblin®) 


} h — a) 
{ p2 sinh Ey + (g + « p2) & cosh & y} cosh & ( 
Cu.) = E{p2cosheh + (g +¢p)& sinh & h} an 
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forO<y<a. Forx« <y<h, yandéareto be inter-changed in the expression 
(3.7). Using the inverse Weber transform formula (3.3) we obtain 


h 


| G (y, a) U(«, p) dx dé. 
0 4 (a8) + Ye (ak) 2 





re) 
e 2 A (r, =) 
G (r, ys p) = = e. 


Expressing the Bessel functions of the first and second kinds in terms of the Hankel 
functions functions and rearranging G we obtain 





: __ 1 f 8,8 coshE(h—y) we 
h — 
x i cosh — (h —- «) U (a, p) da dé jao( 3.8) 
where 
Hy Er) HO Er) 
Bir, §) = Mj ee Os re 
AE@i —  A ER) 
and 
= géEsinh EA 
a eLEDAEy os en ELD 


D (§&) = coshEA + EesinhEh 
E (6 y) = sinh E y + Ee cosh y, 


Taking Laplace’s inversion we obtain 


a h 
Seis ad B(r, &) 
9 (r, y, t) ni [ E D(é&) | | EUG) 
0 0 
t 
hé(h—yp 
a Rd Ahearn uae Al oe y) | U (a, t) sin h(t — +) ds | 
0 
x cosh § (h ~ a) dz dé. pada) 


(3.10) gives the general result for the potential function due to a vertical cylindrical 
wave maker with prescribed time-dependent normal fluid velocity U(y,t). The 
depression of the inertial surface at time ¢ can be obtained from the relation 


1 @ 
C(t) = ez oF & — © ey) (r, 0, 2), scale 
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4. TimMkE-HARMONIC WAVE-MAKER AND STEADY-STATE DEVBLOPMENT 


For a time-harmonic wave-maker we take 
U (y, t) = U(y) sin st (4.1) 


where ¢ is the circular frequency. Then (3.10) gives 


h 


toa T B(r, &) . u cosh £ (h—y) 
9 (r, », t) Ae J vw | oa | zsinot + 


sino t — i t 
i pane fa eso et | cosh & (h — «) d& da. 
...(4.2) 


To find the steady-state development in (4.2) we follow the method used by Rhodes- 
Robinson’. Two cases are required to be investigated according as the integrand of 
the inner integral in the second term of (4.2) hasa pole in the range of integration 
—>0 ornot. Now p2 — o? or equivalently & sinh & h — K* cosh € h has a zero at 
F = &, say, for—& > Oif0 <e« K < | and none if « K > 1 where K = o2/g and 
K* = K(1 —« K)-1._ The latter case is physically interpreted as the inertial surface 
to be “too heavy” while the former as “light”. The two cases are now dealt with 
separately. 


For 0 < e K < 1, we introduce a Cauchy principal value at Ee = Eq (ie. wp = 9) 
and write the inner involving sin p» ¢ in (4.2) as 





— D(&) sinh A p2 — o2 
0 
(g/«)1!2 ; : 
4 cosh &’ (h — y) cosh &' (h -- a) Bir, &') | sin pt ‘a 
= ee P (&') (uw + 9) weth s 
0 
(g}2) 1!” 
a, B(r, £0) cosh £0 (h — y) cosh &0 (= 3)" | LL ALS 
P (&o) ; p— so 


where 
P (é’) = sinh 26’ h + 28’ A and » = # (’). 


By Riemann-Lebesgue lemma the first term iS transient and the integral in the second 


term becomes 7 cos 5 ft as t ~ o9. Thus as ft > ce, we obtain after simplification 


maser ( 3.8 {€ cosh & y — K* sinh & yi 
hE Ft) ee ie ACE 





0 0 
(equation continued on p. 519) 
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x cosh — (h — a) dE U (a) da 


2cosct  B(r, &) cosh Eo (h — y) z Ao 
a seat ee P (Eo) ...(4.3) 


where 
A (&) = & sinh EA — K* cosh Eh } 
h \ ...(4.4) 
Ao = J cosh & (h — a) U (a) du. 
: J 
(4.3) has the alternative representation 


4 Ao cosh Eo (h — y) 
CS P (0) H (& wy) 


Ko (Gn r) cos En (h — y) 
— 4sin ot > An Ki Ki Ena) ay" “slit Beg eae x, A sant. 5) 


{F (§0 r) cos ot + G (Eo r) sin of} 


where &n’s are the solutions of the transcendental equation 


§n sinEn h + K* cosEnh =0, n= Fey ...(4.6) 
and 
A 
An = ! cos En (h — «) U (a) du ...(4.7) 
F (or) = Jo (G0 r) Ji (Eo a) + Yo (8 r) Yi (& a) } 
G (Zor) = Jo (Eo r) ¥1 (Eo a) — Yo (0 r) Ja (Eo a) , (4,8) 
H (Eo a) = Ji (Eo a) + Y? (Ey a) | 


and Ko, K; are modified Bessel’s function. Hence using (3.11), the depression of the 
inertial surface as tf -+ co becomes 





co 
46 Ko (En r) cos En h 
(An) an Tak ac een ve ») As ey Teh 
( g - Ki (§n @) 2énh + sin ZEnh 
Rees AER Ao cosh &o A Fcosst—Gsinor . 
g sinh 2&9 h + 2Egh H (Eq A) 


---(4.9) 
As r > co, this gives 


_ 40 1—e K* 2 1/2 
ea g H(&a) 4 (Er) 


(equation continued on p. 511) 
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= | ¥1 Goa) sin ( Eo — 7 dit ) 


7 
+ I (Eo a) cos ( Egor — rs of )I ...(4.10) 
(4.10) represents outgoing waves at large distance from the wave-maker. 


For « K => 1, there is no pole of the integrand in the second term in (4.2) and 
thus by Riemann-Lebesgue lemma the integral involving sin yt is wholly transient and 
hence f —~ oo 

h ere) 
sin cf B (r, &) {& cosh Ey + ko sinh Ey} 
r,yt)h~ — — |JU ; 
9(r,¥ 1) (+) | E (& sin Eh + ko cosh &A) 





x cosh & (h — «) dE da pl aett) 
where 


re Oke — iy = ...(4.12) 


This has the alternative representation 


= Ko (Snr) cos Sn (h — y) 
g(r, y,t)~ — 4sin st > Ki (Gna) 2%nh+sin2tnh 


n=1 
h 
x J cos Cn (h — a) U (a) da ...(4.13) 
0 
where Cn’s satisfy 
tn sin Cn h — kyo hcos Sn h = 0, n = 1, 2, -.. ...(4.14) 


Then the inertial surface depression as t + oo is 


aoe 40 Se SEE i a Earns 
©O~ FeK—-) © Ki (Gna) 2onh + sin bah 
i! 


h 
x j cos Cn (h — «) U («) da. ...(4.15) 


As r > ©, © (r,t) > 0. Thus when the inertial surface is too heavy a time-harmanic 
disturbance on the wave-maker cannot propagate at large distances from the wave- 
maker. 


5. CONCLUSION 


The problem of a vertical circular cylindrical wave-maker immersed in 4 liquid 
of finite depth with an inertial surface is solved by the use of Laplace transform in 
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time anda suitable Weber transform in the radial co-ordinate. The steady-state 
development of the depression of the inertial surface due to a time-harmonic vertical 
circular cylindrical wave-maker is deduced for a ‘light’ as well as a ‘heavy’ inertial 
surface. In the absence of inertial surface, the results for a time-harmonic wave-maker 
are recovered which can also be duced from Rhodes-Robinson’s’ results in the absence 
of surface tension. 
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The problem of unsteady laminar flow and heat transfer in an incompressible 
micropolar fluid past an infinite porous flat plate has been examined. The 
flat plate is subjected initially to a constant suction velocity followed by a 
step function change for time r’ > 0. An explicit finite difference scheme has 
been used to solve the governing equations of motion and energy. The velo- 
city, microrotation and temperature distribution have been displayed through 
graphs for various values of microrolar parameter R at different time levels. 


1. INTRODUCTION 


Considerable attention has been paid by the researchers on the phenomenon of 
boundary layer flow over a permeable surface through which the fluid is either sucked 
or injected, because of its practical applications to boundary layer control and thermal 
protection in high energy flow by means of mass transfer. A class of unsteady solution 
of Navier-Stokes equations which possesses boundary layer character, is obtained 
when the velocity components are independent of longitudinal co-ordinate with free 
stream velocity either a function of time or constant. 


According to Stuart2, this special consideration leads to the exact solution of 
the flow equations for the arbitrary free stream velocity. 

U(t) = Oot! + f(t)] (1) 
ata very large distance from the solid boundary. Watson® investigated the plane 
flows with special forms of function f (t) in the arbitrary external velocity. Further 
f (t) = 0 in the free stream velocity leads to a simple solution to the plane flows. 
Using this consideration, the problem of unsteady flow with step function change in 
suction velocity have been analysed in Newtonian fluid?. 
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The genera! theory of micropolar fluid which has been formulated and presented 
by Eringen®’6 is deviating from that of Newtonian fluid by accomodating two new 
variables viz. microrotation i.e. spin and micro inertia describing the distributions of 
atoms and molecules inside the fluid elements. Peddieson and McNitt? and Wilson® 
first developed the boundary layer concept for such a fluid and investigated the flow 
past a flat plate. Dey and Nath® studied micropolar fluid flow over a semi-infinite 
plate using parabolic co-ordinates to consider the flow regime including the leading 
edge. Corla ef al.1° solved the steady state heat transfer in a micropolar fluid flow 
over a semi-infinite plate using similarity variables. Chawla!1 obtained the solution 
for unsteady micropolar fluid flow past an infinite plate. 


In the present study, the effect of step function change in suction velocity on 
the flow and heat transfer in an incompressible micropolar fluid past a flat plate has 
been studied by considering the free stream velocity to be constant for all time ¢’ > 0. 
At time t’ = 0, we assume that there isa steady flow over the plate with constant 


suction vs <0. For time t’ > 0, a step function change is made in suction velocity 


which is responsible for unsteadyness of the flow and heat transfer. The velocity and 
temperature 2re assumed to te functions of transverse co-ordinate y’ and time t only. 
An explicit finite difference scheme is employed to solve the governing equations. 


2. MATHEMATICAL FORMULATION 


The unsteady two-dimensional flow of an incompressible micropolar fluid past a 
semi-infinite porous flat plate is considered. The Y’-axis is chosen along the plate 
with leading edge as the origin and Y’-axis, at right angles to it. Let uv’ and v’ be the 
velocity components parallel and normal to the plate respectively. All the fluid pro- 
perties are assumed to be constant. At time t’ = 0, the flow of the fluid is assumed to 


be steady with constant suction velocity ¥ < 0, normal to the plate. At time 1’ > 0, 


the suction velocity is suddenly changed into v, (v < 0 for suction, vs > 0 for in- 


jection) which causes the flow unsteady. At all times, the fluid free stream velocity is 
assumed to be constant and paralled to the plate. Since the resulting flow is super- 


imposed and weak, and the plate is of infinite extent, concerned dependent variables 
may be considered as function of y’ and t’ only. 


The governing equations of momentum and energy for heat conducting micro- 
polar fluid® 12 are 


dy’ 0 
aie ...(2) 
du’ ou’ 02 u' ON’ 
Se ian SAR BO an 
of an’ v ay’ (1 + K) Oy? = K ay’ sasha) 
= 


ey" ...(4) 


UNSTEADY FLOW AND HEAT TRANSFER 5158 








oN’ ae GN 62 N’ ou’ 
ee af at N ae 
U (t) = Ug (constant) ...(6) 
eee er ie a8 K\ ( eu’ \2 
Penge HY ay |= aoa + (H+ 3) (oer) 
KT éu’ _ i oN’ |} 
erry Es + 2N ] 33 lS ] med) 


where P is the density, “ the coefficient of viscosity, cp the specific heat at constant 
pressure, 7 the temperature, y and K are the micropolar material constant and ky the 
thermal conductivity. In view of (1), we take 


— 1, fort’ = 0 


= —}, fort’'> 0 (8) 
where A is the suction parameter. 
Introducing the following non-dimensional variables 

eo" eae E 


u 
= ———;_ t= ————__, *4= 
4 (» + K) (« + K) ties 


ee NwtKk) , ei ieee ON Ter Toc) 
P Uc a | 


where Tw is the temperature at the isothermal wall and 7... the free stream tempera- 
ture. Equations (3), (5) and (7) can be written as 








ie ae oi at 2 8 ON (9) 

ot dy —éOY 1+R oy 

oN oN A @N 2R (i + R) ( du 

aes, ital a a —_ = —I|N — ...(10 

at A ay 1+R ay? R ot oy he 
é 

“2 ig ep al 

ot oy Pr(l+R) 4)? 


AR’ aN \2 du \2 , 2R (2 a1 net) ] 
z E[ ere ( sy ) 2 es a ay 


soe ban 
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where R (= K/p), A (= y/uj) are the micropolar parameters, Pr (= ucp]ks) the 
Prandtl number, E£ (= ue Icy (Tw — Too )) the Eckert number and Re (=Plyv, | j2/?)/ 
is a dimensionless constant. The initial conditions correspond to the solution of the 
steady state problem. 
The boundary conditions can be written as 

u (0, t) = 0, N(0, 1) = 0.0, 6 (0,1) = 1, aty = 9, 


u (co, 1) = 1, N (ce, 1) = 0.0, 6 (09, t) = 0, at y > 0, (12) 


3. MgrHop OF SOLUTION 


The solve the eqns. (9—11) subject to (12), the explicit finite difference scheme 
has been employed. As in Soundalgekar!4, the boundary conditions at y = oo, have 
been taken to be satisfied approximately at y = 6.0. The flow regime defined as a 
semi-infinite strip in time, bounded by y = 0 and y = 6.0 is divided into a grid by 
lines parallel to y and ¢ axes (Fig. 1). The length step Ay is chosen as 0.1 while the 
time step Ar is taken to be 0.0025 so as to satisfy the stability condition as per 
Ralston)». 





jel oe AL OS oe ears ass +——-- 
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j t+-------++-4 ae Re 
! 
nich sag Nn De] 
patie ape Shc eM 
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\ | 
| ae 
\ ! 
I | 
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y—> 


UNSTEADY FLOW AND HEAT TRANSFER 517 


Replacing the time derivatives by forward difference and the space derivatives 
by their corresponding central differences, eqns. (9) — (11) can be re-written as 





yes ae fot | ( Mol tt | Mig — uta) 
(Ay)? 4 2Ay 














epee eed — Nia): 
(1 + R) 2Ay 4 (13) 
= A (Niu, — 2Nt + Nt-1) (Nix. — Ni- 
N* =™M+t ar| Ay ELLE We Oe 2 41 i-1) 
' (1 + R) (Ay)? +A 2Ay 
2R (1 +B) Ui iN 
— UN Shad RE SEA ihe 
R faces MeeL IT )| (4) 
é 
Rie weet 0441 — 261 + Gi_1 i414 — Gi-1 
Oy ets tallest CER (Ay2 fhe 
AR? 
a Fess Nit ey . ( Uig1 — a 
(i + RP Ay 2Ay 
2R 2 ui. — Ut 
+ @+R) (Nels iaee Sayan +) |} (15) 


where 2 <i < 60. 

The initial condition is taken as the steady state solution obtained by solving the 
corresponding steady state finite difference equations through Gauss Seidal method, 
while the boundary conditions (12) transform as 

u, = 0.0, ue = 1.0 
N, = 0.0, Na = 0.0 
6, = 1.0, 621 = 0.0. ...(16) 

The solution of the difference equations is obtained at the intersection of the 

grid lines, called nodes. In Fig. 1, the nodes are specified by double subscripts (i, /) 


with the origin located at the intersection of the lines y = O and? = 0. The value of 
the dependent variables u, N and @ at the nodal points along the lines y = 0 and 


y = 6.0 are known for all time ¢, while unknown values of u¥ (= ui*), Ni (= Ny*?) 


and 9° (= a”), at internal nodes between y = 0 and y = 6.0 for any time t > 9 


are to be determined. 


The non-dimensional shear stress, couple stress and the Nusselt number on the 
wall are respectively given by 
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= eft 2 = (= 
a mi p u | yp’ | oy 1 Ym0 eet lig 
<“* 1 
C My'2' = R2 ( =| 
n= OOOO _ _—————_ = - aie? a! 
(PP tig | vi @yKe + RPA OY Jy .-.(18) 
and 
Nu 00 

(L| vi |p) OY Jym0 (19) 


It has been observed that by increasing the micropolar effects in the fluid, more 
and more heat is transferred from the plate at any time level. 


4. NUMERICAL RESULTS AND DIscussION 
Fixing A = 1.5, Re = 0.5, Pr = 1.0, £ = 0.1 andA = 2.0, the discussion is 
limited to the variation of micropolar parameter R for different time levels viz. =0.1, 
0.5 and 1.0. 


In Fig. 2, ascompared to the Newtonian fluid R — 0.0, the presence of micro- 
polar additives in the fluid increases the velocity throughout. Also as time grows, the 
velocity profiles are further accelerated. The microrotation N is found to have negative 
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Fic. 2. Velocity distribution for different R with f, 
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values for all the parameters under study. Numerically the microrotation increases in 
a region near the plate reaching a maxima with increase in R as well as ¢ (see Fig. 3); 
In Fig. 4, the temperature profiles are drawn. The temperature decreases with in- 
crease in R at all time levels having a square shape profile for R = 5. 
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The difference scheme (13) — (15) is consistent as per condition stated in 
Ralston!5 and Rosenberg!® with the differential equations (9) —(11), which in}turn 
implies that the difference equations actually do approach the governing differential 
equations. The truncation error for the approximation in the velocity, microrotation 
and temperature is [O (At) + O (Ay)?] which tends to become zero as Ay and At tend 
to zero. The programme has also been executed for smaller values of At, viz. At 
— 0,002, 0.001 and 0.0005 and the results compared with those for As = 0.0025, reveal 
no sinfiicant change, thusensuring the convergence of finite difference scheme employed 


here. 
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ELASTIC CYLINDER 
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We consider a random eigen-value problem arising from torsional vibrations 
of acylinder with randomly varying density. The problem is transformed 
into a perturbed Fourier-Bessel equation and the perturbation technique is 
used to obtain its solution. An expression for variance of the eigen-value is 
also derived. 





INTRODUCTION 
We consider the equation 
a2 Ue a2 Ue 1 Ouse ue ue 02 Us 
c { a2 tag Haidar bec ori a = P(r, t) a2 mitt) 


where p is the shear modulus, f the density andr, 6, z are the cylindrical polar co- 
ordinates. This equation arises from torsional vibrations of an elastic cylinder. The 
displacement of the z-axis is assumed to be zero so that 


ug = 0, when r = 0 for all z, ¢. 4 4, 


Also the lateral surfaces of the cylinder are stress free. This gives rise to the boundary 
condition 





or r 
for all z and ¢. . SED), 


This problem has been studied by Achenbach! when the cylinder is hollow and Heath 
and Wood* when the rigidity of the cylinder varies radially. In this paper we use 
perturbation technique to study this problem when the density is a function of r anda 
random parameter +. This problem is of interest due to the fact that many new com- 
posite materials are effectively described as materials with random properties. In an 
earlier paper Zaman® has studied a random eigen-value problem arising from vibration 
of random elastic plates using the variational formulation. However, in this case we 
use perturbation method as the variational method does not seem to be applicable. 


j ; ive itv Pr Pp > 
*Present Address : Department of Mathematical Sciences, King Fahd University of Petroleum 
and Minerals, Dhahran, Saudi Arabia 


Oyo = pb 152 = ue} =OQOatr=a 
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FORMULATION OF THE PROBLEM 
We seek sinusoidal solutions of (1) and put 
ue = V (r) exp [i (wt — « z)] ...(4) 


where w is the angular frequency and « is the reduced wave number. The equation of 
motion (1) and the boundary conditions (2) and (3) transform into 


t . ] 
V(r) + - V(r) + [ P(r, =) @2 — ¢2 — = | V(r) = 0 en &) 


p 
yV(0) =0 
V' (a) = ; V (a) (6) 


where the time dependence ¢‘*! is understood and omitted throughout. Assuming that 
the density P can be written as 


P(r’, t) = (Po + « Pi (r, t)) on) 


with 


where Po is the mean of P and f is a random function with zero mean and € is a small 
parameter. Equation (5) can thus be written as 


rot veortlr-ewm- L]ym=o (8) 
where 
Fe reat BIS 
and 
pout 
rn 


Equation (8) is a perturbation of Bessel’s equation of order]. The boundary condi- 
tions (6) and eqn. (8) form the perturbed random boundary valve problem. Neither the 
perturbed nor the unperturted differential expression is formally self-adjoint and both 
have a singularity at r = 0, 

SOLUTION OF THE PROBLEM 


We transform (8) into Liouville’s normal form by putting 


Vir) = exp | - j f2 dt | u(r). ...(9) 


Equation (8) thus becomes 





u"(r) + [+ - e — ey ©) | u(r) = 0. ...(10) 
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When « = 0, equation (10) reduces to the Fourier-Bessel equation with solutions 


ri/2 Jy (kr), ri/2 Yy (kr), where 


k = (ies a2) 


v- 
which is real and positive whenever r is real and positive (Titchmarsh®; p. 81). 
The transformation (9) transforms the boundary conditions (6) into 
u(o) = Oat r= 0, 
u’ (a) = 2u(a)at r= a. ret) 
The problem (10) with « = 0 is now self-adjoint with one singular end-point at r = 0 


where it is the limit point case®. It has a countable number of real simple eigen-values 
YO, Yl» Y2, --- where yo = 0 and Jj (kn a) = 0 forn = 1, 2, .... 


The associated eigen-functions are 


ug (r) = Ao r3/? 
un (r) = Anril2 Jy (k,r), n > 1. } ed 2) 
With a suitable choice of An, {un}°°, is an orthonormal set in L? (0, a). 
We now write 
Un(r) = un(r) +e Mr) + uM (V) $+... ...(13) 
An = yn + € 7) +e a we .-.(14) 


where un (r) and yn are the nth eigen-function and eigen-value of the unperturbed pro- 


blem and are given by (12). Following Titchmarsh®, we write Tae u as expansions 
in terms of the eigen-functions of the unperturbed problem as 
co 
uD) (r) = >»; anp up (r) 
p=0 
co 
ae (r) a yy Bnp Up (r). Ae R) 


Pp 


tl 


0 


Using (13), (14) and (15) in the perturbed equation (10) and equating the coefficient 
of € to zero, we get 
)” 3 (1) vl) _ vip) 
OE (te aes J + Cre 9 
a) 0. (16) 
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We multiply equation (16) by un (r) and integrate using the orthonormality of un (r)’ s 
to get 


a 
a 2 
rg a foe" | Px (r,t) uw? (r) dr ada 2) 
0 


where we have substituted value of v (r) introduced earlier. In a similar way, we mul- 
tiply (16) by um (r) and integrate to get 


a 





ie 2 
ie ces —— iP | ee Baik 
0 
a 
um (r) dr + | oe un (r) um (r) dr | ...(18) 
0 
i): 
where Y, iS given by (17). 


ESTIMATE ON THE VARIANCE 


We note that yo is a weakly correlated random process (Boyce?). Itis of 


the form 
a 
_() 
i, = \s (r) P (r, t) dr ...(19) 
0 
where 
2 
f() = LBAR i (r) 
pe 
and 


P(r, t) = Pi (r, 7). ...(20) 


The mean square of such a process is given by 


aye aa 
<%q > =| [ Kear f0a) £(r2) dra dre Ql) 
0 Oo 


where K (ry, rz) is the correlation function for the random quantity P1 (r, t). There are 


well documented experimental methods for determination of the correlation function 
appearing in (21) (Corson4 Miller5). Zaman? has given an elementary method to 


derive the correlation function for a statistically homogeneous two phase composite. 
For two phase composites, K (rj, ra) is given by 
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K (ri, r2) = Pi, Pur (A, B) + Past P12 Pre (A,B) 


+ Pi,2 P1,1 Poi (A, B) + pe P22 (A, B) tee) 


where the property 1 is related at two points A and B lying anywhere in the material. 
Pij (A, B) is the probability that point A is in phase i and point B is in the phase j of 
the material and P1,; is the value of P; in the /th material for j = 1, 2. 


Using argu- 
ments based upon elementary probability theory we find’ 
2 —3S | ri — re | 
K =: pS LS Trai TG le 
(ri, r2) <P; > exp | ac Gee) ateo 


where c is the volume concentration of the material forming phase 1, S is its surface 
area and V is the total volume. Thus we can assume the following form 


K (ri, ro) = < py > exp {—# |r — 21} (24) 
where « is a parameter depending upon nature of the composite material. 
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The steady flow of a viscous incompressible conducting fluid beetween two 
coaxial rotating porous cylinders with the outer cylinder bounded by a 
permeable bed is considered. The induced magnetic field is neglected and 
onthe porous bed the boundary condition of Beavers and Joseph is applied. 
An exact solution of the governing equations is found. The velocity are 
temperature distribution of the fluid in free region and in porous region are 
evaluated in dimensionless form. The results are discussed numerically. 


1. INTRODUCTION 


Flow through and past porous media have wider range of applications in many 


branches like chemical engineering, petroleum engineering, soil mechanics and bio- 
medical engineering. 


The study of the flow of a viscous incompressible fluid between two coaxial 
cylinders was first undertaken by Couette? with a view to measuring the viscosity of 
the fluid. Sinha and Choudhary and Jain and Bansal® consider the flow of a viscous 
incompressible fluid betwen two coaxial rotating or non rotating cylinders with the 
walls of the cylinder being either solid or porous. Jain and Bansal® considered the 
flow of a viscous incompressible fluid between two coaxial rotating porous cylinders. 
Jain and Mehta‘ obtained exact solution in a closed form of the hydro magnetic 
equations for an incompressible viscous and electrically conducting fluid flow through 
an annulus with porous walls in the presence of a transverse radial magnetic field. 
Syam Babu’ discussed the flow of a viscous incompressible conducting fluid between 


two coaxial rotating porous cylinders under the influence of a uniform radial magnetic 
field. 


In this paper we study the flow of viscous conducting incompressible fluid bet- 
ween two coaxial rotating porous cylinders with the outer cylinder bounded by a 
permeable bed. The cylinders are rotating with an angular velocities ®1 and ws res- 
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pectively in the same direction. The fluid is injected/sucked at the inner cylinder with 


a constant velocity. The flow in the annulus and in the porous medium is governed 
by the same pressure gradient — dp/dr. 


(u, v, 0) 
(ui, 1, 0) 
(Hr, 0, 0) 
(0, 0, Ez) 
p 

y 

Cp 

7 

To 

he 


2. NOMENCLATURE 


dimensionless fluid velocity in the free region 
dimensionless fluid velocity in the porous region 
magnetic field strength 

electric field strength 

density of the fluid 


kinematic viscosity of the fluid 


= specific heat at constant pressure 
= thermal conductivity of the fluid in zone 1 


= temperature in the fluid region 


ambient temperature 


= heat transfer coefficient 


coefficient of viscosity 


fluid pressure 


= porosity parameter 
= permeability of the porous bed 


= dimensionless constant 


angular velocities of the inner and outer cylinders respectively 
prandtl number 

peclet number 

nusselt number 

magnetic permeability of the fluid 


electrical conductivity of the fluid 


= dimensionless parameter 


magnetic parameter 


dissipation function 


528 K. JAGADEESWARA PILLAI ET AL. 


ri, re = radii of the inner and outer cylinders respectively 
A = suction/injection parameter 
o = slip parameter. 


3. FORMULATION OF THE PROBLEM 


We consider the steady flow of a conducting viscous incompressible fluid between 
two coaxial rotating porous cylinders composed of an insulated material. The cylinders 
terminate at perfect electrodes which are connected through a load. The walls of the 
cylinders being porous witb the outer cylinder bounded bya permeable bed. (Fig. 1) 
The problem is divided into two zones: free zone and porous zone. In zone 1, the 
fluid is governed by the magneto hydrodynamic equations and in zone 2, the flow is 
governed by the Darcy’s law. The solutions in two zones are matched at the interface 





Fic. 1. Physical model. 


between zone | and zone 2 by assuming the slip velocity boundary condition of 
Beavers and Joseph!. A uniform magnetic field Ho is applied in the radial Sasa 
throughout the flow region. The porous zone is assumed as isothermal and so the bed 
is kept at constant temperature 79 which is ambient temperature. The boundary con- 


dition of Rudraiah et al 5. j 
-°. 18 used for temperature at the j 
7 € interface between the porous 


The governing equation | 
Ss for the steady viscous. incom i i i 
eerie ye y : pressible fluid flow in 
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Zone | 
, dau’ y'2 1 dp d2 y' 1 du’ u’ 
dr’ ce Pde sa ( dr * Lar dr’ 7) 
Pea ii ee d2 y’ 1 dv’ y’ 
ae 1G Tee ( dr2 Psd or ax? ) 
+ (0 pe Ho Ez’ — 6 pe? Ho? ¥’) 
Crew) 
dr’ ag 
d(r’ Hr’) _ 
i —— : 
d( Ez’) 
dr’ : 
mal « {( d2T’ 1 aT’ ) 
Oe = ae trae tt 
where 
du’ \2 ue \2 dy’ Ie. i 
s—>[2{()y+( 5) b+ (a r 
Zone 2 
; K { adp’ ; 
vp =~ K[¥ - 29m, | 
pw 
d (r’ uy ) 
dr’ * 
The boundary conditions are 
ue=m,” =r, 7T' = T, atr’ =" 





uv’ 


I 


} 
du ~ ( aw) 
‘ = 
Saeed VK \T | 
dv’ ‘ ) 
vo = vi + re wr, We (ms | a re 
J 


aT’ * ( Z — To). 
fees | 72 
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...B) 


...(4) 


(5) 


...(6) 


sah) 


...(8) 


...(9) 


...(10) 
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The following non-dimensional quantities are used : 









































; ; oe wu. rs v, ra 
r u’ ry 1 = = 
RE Sa - SA SS y “= y » V1 v ’ 
' ' a 2 
My !2 Voil2 Prry ro 5 Bis | 
1] = Ad) Sa = : = Py2 ’ s= rh’ 04 a ? 
r2 we T.. —To 
One 2 Tr = B Act ro uy Prem Uo 
hog es 7 Algom fas == s a ae as em 
2 / 
Pe =A Pr, Lm tl , Ez = “2 , where Ey — —H¢ Ho 
o K (Tz —To) : 2 
' Ge pe” Ho? r 
_ Ar re _ hery —t' 1 
Hr = Fade VR? Nu= : , M2 = FS 
Using the above non-dimensional quantities, eqns. (1) — (9) reduce to : 
Zone | 
du y2 dp d2 4 1 du u 
W dp ig ae as te: ee =e 
aw _ dy 1 dy -(4 M2 Hr? 
dr rar? r dr r2 r2 ) 
M? HI 
+ ae (13) 
d(ru) _ d (rHr) d (Ez) _ 
Se gees oe Ohare pon aD eee ...(14) 
aT _ @T | 1 aT “f(du\?, (u\ 
OOH are aaa pal? Wa) 
dy y \2 
as ( ee ~.) | aC 15) 
Zone 2 
1 d 
“i1=- - a|% <A —2 Vy Qs ] ... (16) 
M? 88 Hr Ez — 209 uy 
V1 = Sawer ray ve 32 M2 Ar: ‘77 ee «astia) 
d (ruy) 
ene (18) 
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The boundary conditions (10) — (11) reduce to 


u=A,v=Q0,, 7 =0 atr = 1 .. (19) 
du ) 
u= Uubdi, Ps —- = (uv1 — 1) | 
dy 


aco 
v= vo + Qe 3, 7 a ze (vol — v4) atr = 53 20) 
| 
T = Tz, ho =Nu Ts. 
dr 


4. VELocity DIsTRIBUTION 


Solving equations (12) — (14) and (16) — (18) with the help of the boundary 
conditions (19) and (20), we get 


u = Al/r 

ot Gre LD) pd err ...(22) 

uy = tad ...(23) 
g 


M2r83 Ez — 2Br Qs 


r2 o2 a M2 52 ...(24) 


y= 


where 
M2 Ez 
B=A[Il + dao], Ci = a5 9, 
Nye = 3A + 8 (AZ +4 (1 + M2 + Aj”) 


The constants C and D are functions of the physical parameters involving in the 
problem. 


5. TEMPERATURE DISTRIBUTION 


Using the velocity of the fluid obtained in zone 1, equation (15) reduces to 


d2T aT Co ‘ 
ss + (1 — Pel ieee = gd 


2 
2p) — Cz3r "2 —Cqr' 
where 
Co = 4L A’, Co = LC2 (Ni — 1)2, C3 = LD? (N2 - 1)? 
Cy = 2 LCD (Ni — 1) (Nz — 1). 


This is a second order ordinary differential equation whose general solution will posses 
a number of singular points. In order to avoid these singularities the solutions are 
obtained for separate cases by using the boundary conditions (19) — (20). 
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(i) No singularities : [Pe A — 2, 2Ni1, 2N2, A, AA — (1 + M?)] 

Co(1 — r-*) Co (1 — 71) C3 (1 — 1) 
2 (2 + Pe) 2N; (2Ni — Pe) 2Ne2 (2N2 — Pe) 


C4 (i — ra) 
A(A — Pe) 


(ii) Singularities; [Pe = — 2, Pe ~ 2N\, 2NeA, A ~ — (1 + M?)] 


i= 


it + Bo (1 — re); ...(26) 


Cor? logr Crit Pe 
= 5 SuhaSNy CINE oe 


Ca (1 — rA) 
A(A + 2) 


(iii) [Pe = — 2 = A, Pe + 2Nj, 2No, M# 1] 


C3 (1 — r2N2) 
2Ne2 (2N2 + 2) 


+ + Bs (1 — r~2); al erp 


Co(1 — r?1) C3 (1 — 242) C' r-2 logr 

en gh el I ee eal. Pi ene # SE ated — Ed aS Pemesy | LS 

2N1 (2Ni +2) | 2Na Ne + 2) 5 + Ba (r-?—1); 
.. (28) 


T 


(iv) [Pe = — 2 = A, Pe + 2No, 2N; = 0, M = 1, Pe ~ 2Ni4] 


T= C'r-*logr _— Cologr 
2 2 


RSA ee net : 
+ 2N2 (2N2 ca 2) + Bs (r 2— 1); 
age) 
(v) [Pe 4A, A = — (1 + M2), 2Ny = 0, Pe # 2No, 2N4] 


_, Col — r-?) 


T= + Ce log r . Cs (1 — r2%2) 


2 (2 + Pe) Pe 2N2 (2No — Pe) 


C4 (1 — rA) 
~ AG =aPe). + Big (r8? 21): .-.(30) 


(vi) [Pe = A = — (1 + M2), Pe # 2Ni, 2No; 2N; = 0, M+ lorA + — 2] 
que Co (1 — r-2) > Cz logr a C3 (1 — r2N2) 
2(1 — M2) 1+ M2 2N2 (2No + 1 + M2) 


C4 logr r-(14M?) 
1+ M2 


(vii) [Pe ~ — 2, 2Ng, A; Pe = 2Ni, 2Ni, 2N, x 0} 


+ + By (r-Q+M?) — 1); ...(31) 


T= Co (1 — r~?) Co r?e log r Ce (te) 
ee ees OF Nee eee + 
2 (2 + Pe) Pe 2N2 (2N2 — Pe) 
(equation continued on p. 533) 
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C4 (1 — r) 


+ oA Geko ee Os ...(32) 


where By ... Bg are known constants which are functions of the physical parameters 
involving in the problem. 


6. DISCUSSION 


The velocity distribution obtained from eqn. (22) and temperature distribution 
obtained from eqns. (25) — (32) in the presence of porous media are evaluated 
numerically and the results are shown graphically. Figures 2-5 gives the velocity 


d=3 
08 € 
as) 
Vv 
an 
0-4 
»d=0 





- 60 


Az3 


| Q, = Q) = 0,8 = 2M = Landa = 0.05. 


Fig. 2. Velocity profiles for Ci = 
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Fig. 3. Velocity profiles for Cy; = 1,0 = 0.05, Q; = Q, = 1,3 =2. 








12 1-4 A ea Et 2:0 
Fic. 4. Velocity profiles for C; = 1,Q, = Q, = 1, M = 3,8 = 2. 


distribution for various parameters involving in the problem and Fig. 6 gives the 
temperature distribution for fixed values of Peclet number and Nusselt number. 
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Fic. 5. Velocity profiles for « = 0.05, Q, = Q, = 1, C, = 9. 


Figures 2-4 give the velocity profiles in the presence of an electric field for 
various values of suction/injection parameter X. When both the cylinders are stationary 
or rotating the velocity increases as the suction/injection parameter increases. The 
velocity is negative for suction and positive for injection. When the porosity para- 
meter ‘o’ increases the velocity increases for suction and decreases for injection for the 
fixed values of the magnetic parameter. When both the cylinders are rotating, the 
velocity increases as magnetic parameter M increases and for A = — | the velocity is 
negative. For small values of M the velocity is negative (M = 1) and for large values 
of M the velocity is positive. To remove the back flow, the magnetic parameter is to 
be increased. For M = 3, the velocity decreases as porosity parameter increases. 


Figure 5 gives the velocity profiles in the absence of electric field. When both 
the cylinders are rotating the velocity decreases with increasing M and increases as 
A increases. For M = |}, and A = —2 the velocity is negative and for all other values 


it is positive. 
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field. 
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Fic. 6. Temperature profiles for C, = 1, M = 1,2, =0,Q,=1. 


Figure 6 gives the temperature distribution against ‘r’ in the presence of electric 
When the inner cylinder is stationary and the outer cylinder is rotating the 


temperature is maximum for A = —3 and minimum for A = —2 andA = 0. The 


temperature increases with the injection parameter. 
increases the temperature decreases. 


wren > 


6. 


When the porosity parameter 
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